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ABSTRACT 

The work reported here was motivated by concern over the use of smooth heat 
flux gages on the otherwise very rough SSME fuel pump turbine blades. To gain 
insights into behavior of such installations, fluid mechanics and heat transfer data 
were collected and are reported for a turbulent boundary layer over a surface with a 
step change from a rough surface to a smooth surface. The first 0.9 m length of the 
flat plate test surface was roughened with 1.27 mm hemispheres in a staggered, 
uniform array spaced 2 base diameters apart. The remaining 1.5 m length was 
smooth. The effect of the alignment of the smooth surface with respect to the rough 
surface was also investigated by conducting experiments with the smooth surface 
aligned with the bases or alternatively with the crests of the roughness elements. 
Stanton number distributions, skin friction distributions, and boundary layer profiles 
of temperature and velocity are reported and are compared to previous data for both 
all-rough and all-smooth wall cases. The experiments show that the step change from 
a rough surface to a smooth surface has a dramatic effect on the convective heat 


transfer. In most cases, the Stanton number and skin friction coefficient drop below 
the smooth- wall correlation immediately downstream of the change in roughness. The 
alignment of the surfaces has only a weak effect of the Stanton number and skin 
friction coefficient just downstream of the interface. 

It is concluded that use of smooth heat flux gages on otherwise rough surfaces 
could cause large errors. It is recommended that heat transfer data collected in this 
manner be used with caution. 
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Chapter I 
INTRODUCTION 


In the design and analysis of many engineering systems, knowledge of the skin 
friction and heat transfer characteristics for turbulent flow are essential. These 
characteristics are usually affected by the shape and condition of solid surfaces. 
Surface roughness plays an important role in many engineering applications such as 
turbine blades, high performance aircraft, heat exchangers, and piping systems. In 
most instances, these surfaces are rough in the aerodynamic sense. 

The objectives of this report are to describe a set of experiments which 
investigates the effects of a step change from a rough to a smooth surface on heat 
transfer in the turbulent boundary layer and to present some results from these 
experiments. The motivation for this work was concern over the results obtained 
when using smooth heat flux gauges on otherwise rough turbine blades. These gauges 
are usually quite smooth. The quality of the data received from this method can 
frequently depend upon the roughness of the surface being tested. If the surface is 
also smooth, the heat flux measured by the gauge will likely be close to the heat flux 
of the surface itself. However, if the surface is rough then the heat flux measured 
with the smooth gauges could be quite different from the actual heat flux on the rough 
surface. The primary application for this research effort involves the similar problem 
of using smooth heat flux gauges for external heat transfer measurements on 
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otherwise rough turbine blades from the Space Shuttle Main Engine (SSME) fuel- 
pump turbine. 

Taylor (1990) measured and classified the surface roughness on turbine blades 
from TF-39 and F-100 aeroengines. Taylor ( 1989) also made roughness measurements 
on the SSME fuel-pump turbine blades. The blades were determined to be very 
rough. The TF-39 and F-100 turbine blades had an average roughness ranging from 
about 1.5 pm to about 10 pm. This was determined to be very rough since the 
boundary layer thickness is on the order of 1 mm. The SSME fuel-pump turbine 
blades had average roughness heights on the order of 15 pm with the boundary layer 
thickness on the order of 0.5 mm. The heat flux gauges used to measure the heat 
transfer are much smoother than surrounding turbine blade surface. How the heat 
transfer is affected by the step change in surface roughness was the primary 
motivation for this research work. 

Flows with step changes in surface roughness are also important in planetary 
boundary layer flows over a land-sea interface and mountain ranges. There is a large 
group of theoretical and experimental work which pertains to step changes in surface 
roughness for these atmospheric flows. 

Background 

Jacobs (1939) performed some of the early experimental work on this subject 
for fully developed channel flows. Antonia and Luxton (1971a, 1972) presented one 
of the most detailed data sets for zero pressure gradient boundary layer flows. They 
performed extensive velocity, turbulence, and skin friction measurements for both 
smooth-to-rough and rough-to-smooth schemes. Their measurements were made using 
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rib-roughened surfaces with the rib crests on the rough surface aligned with the 
smooth surface. Antonia and Luxton (1971b) also studied the effects of a smooth-to- 
rough transition with the bases of the ribs aligned with the smooth section. This gave 
a transition to a rough surface with the roughness elements upstanding. Schofield 
(1975) performed and presented a series of flow measurements for a step change in 
surface roughness with adverse pressure gradients. He also used a nbbed rough 
surface with the crests of the ribs aligned with the smooth surface. Andreopoulous 
and Wood (1982) performed and presented measurements of velocity profiles, 
turbulence quantities, and skin friction distribution for a flow over a smooth plate 
which was roughened in one narrow strip at mid plate using sandpaper. They 
referred to their alignment as "slightly upstanding". 

There have also been theoretical contributions made by Townsend (1965) who 
developed a theory for velocity and temperature profiles development in the internal 
layer which develops after a change in surface condition. Antonia and Wood (1975) 
and Andreopoulous and Wood (1982) each presented prediction methods based on 

numerical solutions of the boundary layer equations. 

As indicated from the previous listed sources, there is a fair amount of 
experimental and theoretical work that has been reported for developments in the 
turbulent boundary layer after a step change in surface roughness. However, there 
exists no heat transfer data for turbulent boundary layer flow with a step change in 
surface roughness. This is not surprising since there is very little good convective 
heat transfer data for rough surfaces in general. 
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Current Research Effort 

The objective of the current research effort was to obtain comprehensive 
convective heat transfer and fluid mechanics data for a step change in surface 
roughness. As stated, this work was performed to gain insight into the use of smooth 
heat flux gauges on the otherwise rough SSME fuel-pump turbine blades. 

In the full scale turbine blade tests, there are two distinct reasonable limits to 
the location of the smooth surface. The first has the smooth surface aligned with the 
base of the roughness elements, or base-aligned, as shown in Figure la. The second 
is the case of the smooth surface aligned with the crests of the roughness elements, 
or crest-aligned, as shown in Figure lb. For this work, both alignments were 
investigated for carefully matched flow conditions. The effects of the alignment of the 
smooth surface with respect to the rough surface at the rough-to-smooth interface 
were investigated to determine if the alignment has a strong effect on the heat 
transfer and fluid dynamics, and if one alignment would be preferred over the other. 

The work consisted of heat transfer and temperature profile measurements as 
well as measurements of velocity profiles, axial turbulence intensity, and skin friction 
coefficient distributions. All measurements were made for both alignments under 
carefully matched flow conditions. Also, the rough-to-smooth data were compared to 
all-rough and all-smooth wall tests made in the same facility. 

Overview 

This report presents the significant heat transfer and fluid dynamics results 
for a step change in surface roughness. Direct comparisons are made between the two 
alignments. The experimental facility and the operation procedures are discussed in 
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Chapter II. A summary of the THTTF operating systems and information concerning 
the Stanton number and skin friction distributions are included. Chapter III presents 
the experimental heat transfer results for the rough-to-smooth surface and compares 
the heat tr ans fer results between the two alignment cases. Likewise, Chapter IV 
presents the experimental fluid mechanics data for the rough-to-smooth surface and 
compares the similarities and difference between the two alignments. Chapter V 
presents the discrete element model and compares the results of the computations 
with the data. A summary and conclusions are presented in Chapter VI. Appendix 
A details the method used to experimentally determine the Stanton number and its 
corresponding uncertainty analysis. Tabular data for the Stanton number and skin 
friction distributions are presented in Appendix B. 
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1.49 m 




d Q = 1.27 mm 
k = 0.635 mm 
L = 2.54 mm 


Figure 1. 


Descriptions of the rough- to-smooth test surface alignment schemes 
and the roughness shape for the h emisph ere roughened surface. 


CHAPTER II 


EXPERIMENTAL FACILITY AND MEASUREMENT TECHNIQUES 

The Turbulent Heat Transfer Test Facility (THTTF) was designed and built to 
study the characteristics of heat transfer and fluid dynamics in turbulent flat plate 
boundary layer flow. The THTTF was originally designed by Norton (1983). It is 
geometrically similar to a test apparatus used at Stanford University to investigate 
turbulent boundary layer flow and heat transfer in flow over a single, porous rough 
surface [Healzer (1974), Pimenta (1975), Coleman (1976), Ligrani (1979)]. 

General Description 

Schematically shown in Figure 2, the THTTF is a closed loop subsonic wind 
tunnel which can deliver a uniform air flow over a set of individually heated flat test 
plates which abut together to form a continuous 2.4 m long test surface. Each plate 
is fitted with its own heater pad and can be maintained for either a uniform 
temperature or uniform heat flux case. The local Stanton number is obtained by 
applying an energy balance to each plate. Boundary layer velocity and turbulence 
profiles are determined with hot-wire anemometry techniques. Skin friction 
coefficients are determined by using hot-wire anemometry or the Preston tube method. 
Thermal boundary layer temperature profiles are measured using a thermocouple 

probe. 
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Descriptions of the apparatus, measurement procedures, and system debugging 
and qualification are discussed in great detail by Coleman et al. (1988) and Hosni et 
al. (1989). The determination of the Stanton number from the measured variables 
and the un certainty associated with the experimentally determined Stanton number 
particular to this work are presented in Appendix A. 

The THTTF is composed of four primary systems: (1) the air flow system, (2) 
the test plate system, (3) the cooling water system, and (4) the automated data 
acquisition and control system. A discussion of each system follows. 

The Air Flow System 

The air flow system is a closed loop design that delivers a uniform velocity, low 
freestream turbulence intensity, controlled temperature air flow at the 10 by 51 cm 
inlet to the 2.4 m long test section. The air velocities can range from 6 m/s to 67 m/s 
which correspond to x-Reynolds numbers from about one million up to ten million at 
the end of the test section. The blower is a Buffalo Forge 45AW industrial blower and 
is driven by a belt and pulley with an Eaton eddy current clutch using a Dynamatic 
18.6 kilowatt electric motor. The air exits the blower and travels through the 1.2 by 
0.6 m overhead duct before entering the header. After turning through the header, 
the air is filtered through a linen particulate filter before passing through a Trane air 
to water 4 row cooling coil. After exiting the cooling coil, the flow is directed through 
a 3.8 cm thick aluminum honeycomb flow straightener and then through a series of 
4 woven screens. The air then enters a three dimensional nozzle which is designed 
to smoothly accelerate the flow without separation [Healzer (1974)]. The flow is 
delivered to the test section at the 10 by 51 cm nozzle outlet with a uniform velocity 
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field and is passed over a 1 mm by 12 mm boundary layer trip immediately upstream 
of the test plates. 

The rectangular cross sectioned test section is fabricated of 1.3 cm thick clear 
plexiglass top and side walls. The bottom wall is the test surface which is comprised 
of the 24 test plates abutted together to form a continuous flat surface. One of the 
sidewalls is equipped with 24 static pressure taps spaced 10.2 cm apart and positioned 
2.5 cm above the center of each test plate. Circular access holes for test probe entry 
are centered over each plate along the center-line of the top wall. These access holes 
are plugged with precision machined stoppers when not in use. The top wall elevation 
is adjustable to maintain the prescribed zero pressure gradient and constant 
freestream velocity in the axial direction. An inclined manometer with a resolution 
of 0.06 mm of water is connected to the sidewall static pressure taps to measure the 
pressure gradient in the flow direction. For the rough to smooth tests the static 
pressure variations along the test section were less than 4.2 mm for a freestream 
velocity of 66 m/s and 0.19 mm for velocity of 12 m/s, with corresponding maximum 
pressure coefficients of 2.0 x 10' 5 and 2.6 x 10‘ 5 respectively. 

After exiting the test section, the air passes through an adjustable plexiglass 
diffuser and a series of screen inserts and enters a wooden vaned diffuser before 
entering the blower plenum. Filtered make-up air is added in the blower plenum to 
replace air lost through leakage in the air conduits. 

The blower and motor are affixed to a massive concrete pad with vibration 
damping feet to minimize vibrations. The test section is isolated from key vibration 
sources through the use of flexible couplings at key locations. The noise in the 
overhead duct, plenum, and header is damped through the use of batt insulation 
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covered by rigid fiberglass insulation board in these components. Also, the blower and 
motor are housed inside vented, insulation lined boxes for noise abatement. 

The Plate System 

The bottom surface of the THTTF test section is composed of individually 
machined plates that fit together to form a continuous flat surface. As shown in 
Figure 1, the first 0.91 m of the test section was composed of 9, 10.2 cm wide plates 
roughened with 1.27 mm diameter hemispheres, and the remaining 1.49 m length was 
made up of smooth plates. The roughness elements were spaced 2 base diameters 
apart in a staggered array as shown in the figure. Each rough plate was machined 
from a solid aluminum blank, so the contact between the roughness element and the 
smooth portion is perfect. The portion of the test surface immediately downstream of 
the rough-to-smooth interface was composed of 4 smooth plates with 2.54 cm width. 
These were installed to give better resolution of the Stanton number behavior after 
the step. The remaining smooth-wall plates were 10.2 cm wide. Dowel pins, as shown 
in Figure 3, are used to secure the plates together. The allowable mismatch at the 
joint between two plates is 0.013 mm. The plates rest on support rails that are heated 
to the plate temperature to minimize conduction losses. The test section support 
structure is fully encased in a removable wooden shroud which is amply insulated to 
prevent heat transfer with the laboratory environment. Figure 4 shows a cross section 
of the THTTF test section. 

The hemispherically roughened plates with L/d 0 = 2 were chosen because on 
the 2.4 m test section plate, the roughness scales reasonably well with the roughness 
on the SSME fuel pump turbine blades. Taylor (1989) made comparisons of the 
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THTTF test section to the SSME blades. Assuming the heat flux gauges to be 
1.25 mm in diameter, the following table compares the conditions in THTTF test 
section to estimated conditions in ground tests on the actual SSME blades at Calspan 
in Buffalo, NY. 



IV8 

k/8 

St/St* 

Re x 

SSME 

4 

0.003-0.03 

= 2 

= 10 6 

THTTF 

1-2 

0.025 

= 2 

= 10 6 

L is the length 

of the smooth heat flux gauge on 

the SSME 

turbine blade (1.25 mm) 


or the length of the four combined 2.54 cm wide THTTF plates (10.2 cm), k is the 
roughness height, 5 is the boundary layer thickness, and St s is the smooth surface 
Stanton number. The SSME Reynolds number is based on estimated test conditions 
at the mid chord. There is good agreement between the scaling of the THTTF test 
section and the actual SSME blades. 

Each plate in the THTTF contains two thermistors for temperature 
measurement and is fitted with a flexible, rubber-encased, electric heating pad. The 
plate temperature is measured to within ±0.1 C using the two thermistors which are 
mounted in wells drilled into the lower surface of the plate. The power to each plate 
is regulated by a motorized transformer, which is controlled by the Automated Data 
Acquisition and Control System. Experience in performing heat transfer tests in the 
THTTF has shown that plate temperatures can be held to within ±0.1 C of a desired 
constant wall temperature for the entire 2.4 m length of the test surface. Design 
computations have shown that under these heating conditions the plate and the 
associated roughness elements can be considered to be at a uniform temperature. 
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The Cooling Water System 

The cooling water system is needed by the THTTF to remove excess heat from 
the air and maintain the air at some constant temperature. This system is composed 
of a 570 liter water storage tank, a Trane 4 row air-to-water cooling coil, a motorized 
ball valve, a centrifugal pump, and PVC piping to circulate the cooling water. The 
cooling water is pumped from the storage tank to the cooling coil. Heat is transferred 
to the cooling water from the warm test air. A proportional amount of the warmer 
water may be dumped into a floor drain through the motorized ball valve. The ball 
valve is manually adjusted so that the water temperature in the loop and, thus, the 
air temperature in the test section is held at some constant value. The water level 
in the storage tank is held at a desired level with make-up cooling water from the 
building supply by depth sensing valves located in the storage tank. 

The Data Acquisition and Control System 

The THTTF is monitored and controlled by a Hewlett-Packard Series 9000 
Model-220 personal computer and a Hewlett-Packard 3045A Automatic Data 
Acquisition and Control System (ADACS). This includes an HP-3437A high speed 
system voltmeter, an HP-3456A high resolution digital voltmeter, and HP-3497A data 
acquisition/control unit and a number of special function plug-m assemblies. 
Suryanarayana (1986) gives a detailed discussion of the ADACS and its use in the 

THTTF. 

Information is relayed to the personal computer via transducers that are wired 
into ADACS and are used to monitor the THTTF. Using the operating condition 
information, the personal computer decides on the proper response based on 
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programming, and sends the properly controlled commands to the ADACS. The 
ADACS can control the rail heaters, the plate heaters, the cooling water dump valve, 
and the blower motor. When the proper equilibrium conditions are met, the personal 
computer instructs the ADACS to collect the necessary data. 

Measurement Techniques 

The techniques used to determine the Stanton number, the mean velocity and 
turbulence quantities, the skin friction coefficient distributions, and the thermal 
boundary layer temperature profiles are discussed below. 

Stanton Number Determination 

By applying an energy balance to each test plate in the test section, the data 
reduction equation for the experimentally determined Stanton number is found to be 

st - (i) 

p U„AC p (T w -T 0 ) 

The plate heater power, W, is measured using a precision wattmeter. Using a gray- 
body enclosure model, the radiation heat loss, q r , is determined with the emissivity 
of the nickel plated plates estimated to be 0.11. The conductive heat loss rate, q^, is 
calculated using an experimentally determined value of overall conductance from plate 
to rail, (UA) pr , and an experimentally determined value of overall conductance from 
plate to plate, (UA) . The conduction losses are minimized by insulating underneath 
the test plates and actively heating the support rads. The density and specific heat 
of the air are determined from moist air property data using measurements of 
barometric pressure and wet and dry bulb temperatures in the THTTF test section. 
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The freestream velocity is determined with a Pitot tube and specially calibrated 
pressure transducers. Specially calibrated thermistors are used to measure the air 
and plate temperatures. The freestream total temperature is calculated using a 
recovery factor of R = 0.86 (Eckert and Goldstein, 1976) for the freestream thermistor 
probe. The derivation of the Stanton number data reduction equation from an energy 
balance on a test plate and the details of determination of each of the variables used 
are presented in Appendix A. 

The uncertainty analysis of the Stanton number is based on the ANSI/ASME 
Standard on Measurement Uncertainty (1986), following the procedures of Coleman 
and Steele (1989). The bias limits for the sixteen variables involved in the 
determination of the Stanton number were estimated by Hosni (1989) and are 
presented in Appendix A. Some of the elemental contributions to the bias limits were 
correlated since the thermistors were calibrated against the same standard. These 
correlated biases, which reduce the overall uncertainty in the experimentally 
determined Stanton number, were accounted for in the uncertainty analysis. 

During the design and construction phase of the THTTF, a detailed uncertainty 
analysis of the determination of Stanton number was made. This analysis showed 
that relative to bias limits corresponding to the measured variables, the precision 
limits were negligible. Therefore, replications of the Stanton number at a given set 
point should show negligible scatter. In previous work in the THTTF, it has been 
shown that for low freestream velocities ( ^ 12 m/s for the smooth plates and 

< 6 m/s for the rough plates ) the time constant of the THTTF is large enough 
such that variations in line voltage and cooling water temperature reduce the ability 
to maintain a strict steady state. As discussed by Hosni et al. (1989), a 95 percent 
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confidence estimate of the precision limit in the 6 m/s rough plate and 12 m/s smooth 
plate Stanton numbers was estimated as 3 percent. The 3 percent precision limit 
contribution was combined with the bias limit by a root-sum-square to obtain the 
overall uncertainty in Stanton numbers. 

Mean Velocity and Turbulence Intensity 

Local mean velocities are measured with the horizontal hot-wire. The 
horizontal hot-wire probe is mounted on a micrometer arrangement which allows 
proper positioning along the x, y, and z axes. Using a fourth order least squares 
calibration equation, velocities are calculated from 1000 instantaneous anemometer 
output voltages measured 0.01 seconds apart. For each velocity profile, the probe was 
positioned in the freestream and the reference velocity, U^, was measured. The probe 
was then lowered to a height above the smooth wall as determined by a fixed keel. 
The keel height was set to prevent contact between the hot-wire probe and the plate 
surface. Measurements were taken at intervals such that the normalized velocity 
(u/UJ increased 1-2 percent from point to point in the near-wall region and 2-4 
percent in the outer region. The overall uncertainty in the values of u were estimated 
by Coleman et al. (1988) to be ± 2 percent. 

The measurements of the fluctuating longitudinal velocity component u ' 2 were 
made with the horizontal hot-wire in conjunction with the mean velocity 
measurements described above. The longitudinal velocity fluctuation u ' 2 was taken 
as the variance of the 1000 computed velocities. The overall uncertainty associated 


with the measurements of u ' 2 is ± 5 percent. 
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Skin Friction Coefficient Determination 

The skin friction coefficients over the smooth-wall surface were determined 
using the method of Preston (1954). Preston’s method uses a simple Preston tube 
resting on the wall surface and depends upon the assumption of a universal inner law 
(law of the wall) common to smooth wall turbulent boundary layer flows. The 
difference between the total pressure at the Preston tube of 1.6 mm inside diameter 
and the undisturbed static pressure measured at a pressure tap in the sidewall at the 
same x-location as the Preston tube was measured with the same calibrated pressure 
transducers used to determine the freestream velocity. This difference in pressure 
was used with Patel’s (1965) calibration equation to solve for the skin friction 

coefficient. 

The Preston tube method of determining the skin friction coefficient is valid 
only for smooth-wall turbulent boundary layer flows because calibrations have only 
been made for smooth walls. Unfortunately, the region just downstream of the step 
change in surface roughness does not have a well developed smooth-wall inner layer. 
Therefore, an additional uncertainty for the Preston tube measurements of skin 
friction coefficient made in this region exists. However, as shown later, the inner 
layer develops quickly, and all of the Preston tube measurements, with the exception 
of those made at the first smooth plate, are considered to be valid determinations with 
uncertainties of ±6 percent (Coleman et al.,1988). 

Temperature Profile Measurements 

A Type E (chromel-constantan) butt-welded, 0.25 mm diameter thermocouple 
probe was used to measure the time mean temperature distribution across the 
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boundary layer. The probe was mounted on the same micrometer arrangement used 
for the velocity profile measurements, and the method of traversing the thermal 
boundary layer was the same as the procedure for the hot-wire probe discussed 
previously. The probe was positioned in the freestream and the reference freestream 
temperature was measured. The probe was then lowered to a position 1.19 mm above 
the wall for the initial measurement. Measurements were taken throughout the 
thermal boundary layer until the freestream was reached. The probe output, in 
millivolts, was converted to temperature using Hewlett-Packard software. The 
uncertainty associated with the temperature measurements was determined to be 
± 0.09 °C (Coleman et al., 1988). 
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Figure 2. Schematic of the Turbulent Heat Transfer Test Facility (THTTF). 
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CHAPTER III 


HEAT TRANSFER RESULTS 

This chapter presents the experimental heat transfer results for zero pressure 
gradient, constant wall temperature, incompressible turbulent boundary layer flow 
with freestream velocities of 6, 12, 27, 43, 58, and 66 m/s and a nominal wall 
temperature of 44 C for both the base-aligned and crest-aligned cases. The main 
emphasis of the results is the Stanton number measurements in the Turbulent Heat 
Transfer Test Facility (THTTF) for the two alignments. Furthermore, the behavior 
of thermal boundary layer thickness and enthalpy thickness is also investigated. 

The Stanton number data were determined from energy balances on each test 
plate, as outlined in Appendix A. Using the ANSI/ASME Standard on Measurement 
Uncertainty (1986) and following the procedures of Coleman and Steele (1989), the 
uncertainty in the experimentally determined Stanton number was estimated. As 
discussed in detail in Appendix A, the overall uncertainty for the Stanton number 
data ranged from about ±2 percent to ±5 percent for the 0.1 m plates, depending upon 
freestream velocity. The overall uncertainty in Stanton number associated with the 
narrower 0.025 m smooth plates at the rough-to-smooth interface ranged from ±4 
percent to ±11 percent, depending upon freestream velocity. The plots of Stanton 
numbers presented in this chapter include uncertainty bars on various data points 
which show the typical overall uncertainty limits on Stanton numbers. 
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Temperature profiles of the thermal boundary layer were measured with the 
thermocouple probe, as described in Chapter II, for freestream velocities of 6, 12, 27, 
43, 58, and 66 m/s. The profiles were taken at a location 0.94 m downstream of the 
test surface leading edge. This location is just downstream of the rough-to-smooth 
interface. For the base-aligned tests, additional temperature profiles for freestream 
velocities of 12 and 58 m/s were measured at locations 0.94 m, 1.05 m, 1.15 m, 1.35 m, 
and 1.55 m downstream from the leading edge. These temperature profiles were used 
to deter min e the thermal boundary layer thicknesses. The enthalpy thickness 
corresponding to each Stanton number was determined by numerically integrating the 
energy integral equation [Kays and Crawford (1980)]. 

X 

A 2 = J St(x)dx ^ 

o 

Base-Aligned Case 

Figures 5 through 10 present the base-aligned Stanton number data plotted 
versus Re x for nominal freestream velocities of 6, 12, 27, 43, 58, and 66 m/s, 
respectively. The corresponding x-Reynolds numbers range from 100,000 to 10 
million. The x-Reynolds number immediately downstream of the step ranges from 
300,000 to 3,000,000. Also shown on these plots is the corresponding THTTF heat 
transfer data for a smooth surface over the entire 2.4 m length of the test section 
(Coleman et al., 1988). For the 6 m/s case shown in Figure 5, a line representing the 
smooth wall case for turbulent boundary layers (Coleman et al.,1988) 

St = 0.185 [log 10 (i?e x )]- 2 - 584 Pr-°- 4 (3) 

was plotted for comparison since no smooth-wall data were taken for 6 m/s. These 
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figures show that the step change from a rough surface to a smooth surface has an 
immediate and large effect on the heat transfer. For each of the frees tream velocities 
the Stanton number rapidly approaches the smooth wall case. In most of the cases, 
the Stanton number immediately after the step change in roughness falls below the 
corresponding smooth-wall values and rises back to the smooth-wall values further 
downstream of the interface. 

The U„ = 6 m/s case shown in Figure 5 behaves differently. In spite of the trip 
at the nozzle exit and the rough surface, the flow remains laminar for a considerable 
length. The flow becomes fully turbulent at a Reynolds number of about 200,000, and 
a transitionally rough boundary layer is established for a short distance before the 
rough-to-smooth interface. According to Hosni et al. (1989) this case is transitionally 
rough 1 in the aerodynamic sense while the higher velocities are fully rough boundary 
layers over the rough portion of the test surface. In the region between the rough-wall 
boundary layer and the developing smooth-wall boundary layer, the Stanton number 
decreases rapidly in a smooth, continuous fashion to the new smooth-wall condition. 
For a freestream velocity of 12 m/s shown in Figure 6, the Stanton number drops 
immediately to the smooth-wall condition. There appears to be a slight dip in the 
rough-to-smooth Stanton number below the all smooth-wall Stanton number data. 
Freestream velocities of 27 and 43 m/s shown in Figures 7 and 8, respectively, have 


1. Turbulent flows which are influenced by surface roughness are usually divided 
into 3 regimes. Aerodynamically smooth flows are those where the roughness effects 
are so small that the flow behaves as if the wall were smooth. Fully rough flows ^are 
those where the roughness so dominates the momentum transport to the wall t a 
viscous effects are negligible. In turbulent pipe flow fully rough flows are those where 
the friction factor is no longer a function of the Reynolds number. Transitionally 
rough flows are those at Reynolds numbers between, where viscous and roughness 
effects are both important. 
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an immediate drop in the Stanton number after the rough-to-smooth interface. Both 
cases show a definite dip in the Stanton number below the equivalent smooth-wall 
values after the step change in surface roughness. At the higher velocities of 58 and 
66 m/s shown in Figures 9 and 10, the change in the Stanton number to the smooth- 
wall case is still immediate; however, the dip below the equivalent smooth-wall data 
is weaker than for the 27 and 43 m/s cases. 

Overall, the transitionally rough flow has a smooth continuous change in the 
Stanton number from the rough surface region to the smooth surface region. Fully 
rough cases, however, have an immediate abrupt change in the Stanton number after 
the rough-to-smooth interface, with the Stanton numbers after the interface falling 
below the equivalent all smooth-wall values. As the freestream velocity increases, the 
change becomes slightly more abrupt. 

Figures 11 and 12 show plots of the non-dimensional temperature profiles 
versus y/A for freestream velocities of 12 and 58 m/s for the base-aligned case. These 
plots show that the near wall region of the thermal boundary layer downstream of the 
rough-to-smooth interface rapidly takes on the smooth-wall profile characteristics. 
The outer regions of the profiles further downstream of the interface gradually 
approach the all smooth-wall case. 

Figure 13 shows the temperature data for U,,,, = 12 m/s plotted in T versus y 
coordinates. Also shown on this figure is the law of the wall for the thermal boundary 
layer for a smooth wall as given by Kays and Crawford (1980) 

T* = 2.195 ln(y + ) + 13.2Pr - 5.66 (4) 

The first smooth-wall T + profile immediately after the roughness interface is much 
like the last rough-wall profile. The profiles further downstream of the interface 
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approach the smooth-wall T+ levels. The uncertainties in the T + values are relatively 
large. 

Figures 14 and 15 show plots of enthalpy thickness versus x for the base- 
aligned cases at 12 and 58 m/s. The corresponding all-rough data reported by Hosni 
(1989) are also shown. The enthalpy thickness after the step continues to grow; 
however, the growth rate is less than that a completely rough test surface. 

Alignment Comparison 

A direct comparison of the Stanton number data for the two alignments is 
shown in Figures 16 through 21. For the 6 and 12 m/s cases, results are 
indistinguishable from each other within the data uncertainty. The data at the rough- 
to-smooth interface exhibit the same behavior for both the base-aligned and crest- 
aligned cases. For the 27 m/s case, the Stanton number for the crest-aligned case on 
the first 0.025 m smooth test plate is slightly larger than its base-aligned counterpart; 
however, there is a large overlap of the uncertainty intervals. The same behavior is 
seen at the higher freestream velocities of 43, 58, and 66 m/s, with the first crest- 
aligned Stanton number after the interface slightly larger than its base-aligned 
equivalent with a small amount of overlap of the uncertainty bands. This effect is not 
seen for any test plate Stanton numbers further downstream of the first smooth-wall 
plate regardless of freestream velocity. For the conditions of these experiments, any 
effect of surface alignment is very small and observed only in the region immediately 
downstream of the step change in surface roughness. 

As stated earlier, for the crest-aligned case mean temperature profiles were 
measured at an x-location 0.04 m downstream of the rough-to-smooth interface for all 


26 


freestream velocities. Figure 22 shows a comparison of the mean temperature profile 
at this location for each of the two alignments at a freestream velocity of 12 m/s. 
There is no real discernable difference between the two profiles of either alignment. 
Figure 23 shows the same plot for a freestream velocity of 58 m/s. 
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Figure 5. Stanton number data versus Re x for the base-aligned case for U, 
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Figure 9. Stanton number data versus Re x for the base-aligned case for U M = 58 m/s. 
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Figure 11. Nondimensional temperature profiles versus y/A for U M = 12 m/s for the base- 
aligned case. 
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Figure 12. Nondimensional temperature profiles versus y/A for U M = 58 m/s for the base 
aligned case. 
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Figure 13. Temperature profiles at U M = 12 m/s for the base-aligned case compared with 
the thermal law of the wall. 
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Figure 14. Enthalpy thickness versus x for the base-aligned case for U M - 12 m/s 
compared with the all- rough case. 
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Figure 15. Enthalpy thickness versus x for the base-aligned case for U M - 58 m/s 
compared with the all-rough case. 
















600 ' 



Figure 20. Comparison of Stanton number data versus Re x for the base-aligned and crest 
aligned cases for = 58 m/s. 
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Figure 22. Comparison of the nondimensional temperature profiles versus y/A for the 
base-aligned and crest-aligned cases for U M = 12 m/s. 



45 



Figure 23. Comparison of the nondimensional temperature profiles versus y/A for the 
base-aligned and crest-aligned cases for U M = 58 m/s. 



CHAPTER IV 


FLUID MECHANICS RESULTS 

This chapter presents the fluid mechanics results obtained in the boundary 
layer over the rough to smooth surface for both the base-aligned and the crest-aligned 
cases. All of the fluid dynamics results are for zero pressure gradient, isothermal, 
incompressible boundary layer flows. Profiles of mean velocity and the turbulence 

intensity quantity ~uP‘ were measured with a horizontal hot-wire using the techniques 

outlined in detail by Hosni et al. (1989). These profiles were used to obtain the 
boundary layer thickness, displacement thickness, and momentum thickness along the 
rough-to-smooth surface. The skin friction coefficient distribution was obtained using 
the Preston tube method. All of the data obtained are compared with previously 
reported THTTF data for all-rough and all-smooth surfaces. 

Base-Aligned Case 

Boundary layer velocity profiles and boundary layer integral parameters such 
as the boundary layer thickness (5), the displacement thickness (8^, and the 
momentum thic kn ess (82) were obtained for the rough to smooth surface at nominal 
freestream velocities of 12 and 58 m/s. These profiles were measured at locations of 
0.95 m, 1.05 m, 1.15 m, 1.35 m, 1.55 m, 1.75 m, and 1.95 m downstream of the test 
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section leading edge, for the base aligned case, with the rough-to-smooth interface 
located at 0.90 m. 

Figure 24 shows a plot of the rough-to-smooth velocity profiles in u/U.. versus 
y/S coordinates taken with a hot-wire at a nominal freestream velocity of 12 m/s. Also 
included on the plot is the velocity profile for the last rough plate, x = 0.85 m. This 
plot shows the profile immediately downstream of the rough-to-smooth interface 
quickly assumes the smooth-wall characteristics in the near-wall region, but resembles 
the rough-wall profile in the wake region. Further downstream of the step, the 
profiles gradually assume a smooth-wall shape. Figure 25 shows a composite plot of 
the boundary layer mean velocity profiles taken with a hot-wire at a nominal 
freestream velocity of 12 m/s plotted in y versus u/LL coordinates. The u/U^ abscissa 
is plotted with a multiple origin to show the progression of the velocity profiles 
downstream of the rough-to-smooth interface. The plot also shows the corresponding 
rough-wall velocity profiles under carefully matched flow conditions for the identical 
locations. This figure shows that after the interface, the velocity profile is quick to 
deviate from the rough-wall profile near the wall, however, there is quite a distance 

before the fully smooth wall profile is obtained. 

Profiles of axial turbulence intensity were also determined with the horizontal 
hot-wire at the same locations downstream of the rough to smooth interface as the 
mean velocity profiles. Figure 26 shows the rough-to-smooth profiles of axial 
turbulence intensity normalized by U„ and plotted against the y-position normalized 
by the boundary layer thickness 5. For the 12 m/s profiles, the sharp near-wall peak 
typical for smooth-wall profiles is seen for all x locations except the x = 0.95 m 
location immediately after the step. Figure 27 shows the turbulence intensity profiles 
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for a freestream velocity of 12 m/s plotted in coordinates y versus y u /2 /U„ on a 
multiple origin plot. Also shown on this plot are the rough-wall profiles for the same 
x-locations. The near-wall regions for the rough-to-smooth case Quickly deviate from 
the all-rough scheme. However, the flow still has the rough-wall characteristics 
further out in the wake region for a considerable distance downstream. 

Figure 28 shows a plot of 5 versus x for a freestream velocity of 12 m/s. The 
boundary layer thickness, 5, is taken as the distance above the plate at which the 
boundary layer velocity was within 1 percent of the freestream velocity. This figure 
contrasts the boundary layer thicknesses for the all rough-wall and the rough-to- 
smooth cases. After the change in surface roughness, the boundary layer thickness 
grows at a slower rate than in the all-rough case. 

The boundary layer displacement thickness was determined from numerical 
integration of the incompressible flow displacement thickness definition using the 
mean velocity profile data obtained with the horizontal hot-wire, 

8,-j" <*> 

Likewise, the boundary layer momentum thickness was found by numerically 
integrating the incompressible flow momentum thickness definition using mean 
velocity profile data 



Figure 29 contrasts the boundary layer displacement thickness and momentum 
thickness for the all-rough and rough-to-smooth cases. While the all-rough 
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displacement increased rapidly with x, the displacement thickness for the rough-to- 
smooth case remains more or less constant after the interface before slowly growing. 
After the rough-to-smooth interface, the momentum thickness closely follows the 
equivalent all-rough case for a substantial distance (0.4 m). After this distance, the 
rough-to-smooth case departs from the all-rough case with a slower growth rate. This 
effect is clearly seen in Figure 30. Figure 30 shows the nondimensional shape factor, 
H = Sj/Sg, versus x. There is a noticeable drop in the shape factor after the flow 
crosses the rough-to-smooth interface. The all-rough case has a more or less constant 
value for the shape factor equal to 1.6. The rough-to-smooth case, however, gradually 
falls toward the value of 1.3, the widely accepted value for the shape factor for a flat 
plate in turbulent flow. Antonia and Luxton (1972) reported the same behavior with 
shape factor in their experiments. 

Antonia and Luxton (1971b) discussed, for a step change in surface roughness, 
the existence of two layers in the flow field downstream of the roughness interface. 
The boundary layer, 5, was defined previously to be the distance above the plate at 
which the boundary layer velocity was within 1 percent of the freestream velocity. An 
internal layer, 5j, grows immediately downstream of the step due to the rapidly 
varying shear stress distribution. Shown in Figure 31, the flow outside the internal 
layer remains mostly unaffected by the change in surface roughness, while the flow 
inside the internal layer assumes a smooth-wall characteristic. Antonia and Luxton 
plotted their mean velocity profiles in u/U„ versus y 1/2 coordinates and showed that 
a " kink " or "knee" occurred in the plot at the edge of the internal layer. Figure 32 
shows the mean velocity profiles for 12 m/s plotted in these coordinates with a shifted 
origin. Lines are drawn on the profiles to show quantitatively the location of the 
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"kink" or "knee" point and, therefore, the internal layer thickness. Figure 33 shows 
a plot of the boundary layer thickness, 5, and the internal layer thickness, Sj. The 
internal layer grows more rapidly than the bounda r y layer, since ultimately it will 
envelop the entire boundary layer. At x = 1.9 m the ratio of 5/8 = 0.8. The mean 
velocity profiles shown in Figure 24 are repeated in Figure 34 with the boundary layer 
thickness, 5, and the internal layer thickness, 8j, shown on the plot. There is good 
agreement between the point where the rough-wall and rough-to-smooth wall profiles 
meet, and the location of the inner layer edge found using Antonia and Luxton’s 
method for the first few profiles. For the latter profiles the agreement is not as good. 
This is due, partially, to difficulty in determining the location of the "kink" using 
Antonia and Luxton’s method on these latter velocity profiles. 

In addition to the mean velocity profiles, local skin friction coefficients were 
determined along the smooth-wall portion of the test surface using Preston tube 
measurements. Figure 35 shows the skin friction coefficient distribution for the 
rough-to-smooth surface at a freestream velocity of 12 m/s. Also shown on the plot are 
the all-rough data from Hosni (1989) and the all-smooth wall data from Brown (1988). 
Brown’s data were taken with the same Preston tube. Hosni’s rough-wall data have 
an uncer tain ty of ± 10 percent. The values of C f determined from the Preston tube 
method have an uncertainty of ± 6 percent except for the value at the first smooth 
plate, as discussed in Chapter II. As with the Stanton number, the values of C f fall 
below the smooth-wall values immediately after the rough-to-smooth interface and 
then rise back to the smooth-wall values. Antonia and Luxton (1972) and 
Andreopoulous and Wood (1982) reported a similar trend in their experiments. The 
skin friction coefficient data shown in Figure 35 are tabulated and shown in Appendix 
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B. Data at 58 m/s give similar results and are also tabulated in Appendix B. 

The mean velocity profiles are plotted in the inner variable coordinates, u + 
versus y + , for the freestream velocity of 12 m/s in Figure 36. Also shown is the 
reference "law of the wall" for smooth surfaces [Kays and Crawford (1980)]. 

(7) 

u* = 2.44 lny + + 5.0 

In these coordinates, the mean velocity profiles show how quickly the flow assumes 
the smooth wall profiles. The profile at the x = 0.95 m location displays a slight 
nonlinearity in the log region. The later profiles, however, show an excellent 
agreement with the smooth wall log region. In contrast to the thermal boundary layer 
profiles shown in Figure 13, the near wall velocity profiles almost immediately fall on 
the smooth- wafl line. Recall that the temperature profile 0.04 m downstream of the 
interface almost matched the rough- wafl profile. Since the profiles downstream of the 
rough-to-smooth interface rapidly match the smooth wall log law, using the Preston 
tube method with Patel’s calibration to determine the skin friction coefficient was 
assumed valid as discussed in. Chapter II. 

Alignment Comparison 

Figure 37 compares the mean velocity profiles for the base-aligned and crest- 
aligned cases for a freestream velocity of 12 m/s. Both profiles were measured at a 
location 0.05 m downstream of the rough-to-smooth interface. As seen in this plot, 
there is virtually no difference between the two velocity profiles. There is a slight 
difference between the data for the two cases, but this difference is well within the 
data uncertainty. For the axial turbulence intensity profiles shown in Figure 38, the 
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data are virtually indistinguishable between the two alignment cases. A similar 
comparison for skin friction coefficient, C f , is shown in Figures 39 and 40 for 
freestream velocities of 12 and 58 m/s, respectively. The crest-aligned values of C f 
appear to be slightly less than the corresponding base-aligned data. However, since 
the skin friction coefficient data have an uncertainty of ±6 percent, no real difference 
in the value of C f for the two alignments can be concluded. For the above 
comparisons, any differences in the flow field due to alignment effects are not 
substantial. 
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Figure 26. Axial turbulence intensity profiles for the rough-to-smooth test surface at 
increasing x-locations for U,^ = 12 m/s. 
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Figure 27. Axial turbulence intensity profiles for the rough-to-smooth and all-rough test 
surfaces at increasing x-locations plotted on a shifted axis for U M = 12 m/s. 






Figure 28. Comparison of boundary layer thickness versus x for the rough-to-smooth 
surface and all-rough surface for U,,,, = 12 m/s. 
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Figure 30. Comparison of shape factor versus x for the rough-to-smooth surface and the 
all-rough surface for = 12 m/s. 
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Figure 31. Internal layer and boundary layer regions of flow downstream of a step change 
in surface roughness. 
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Figure 32. Mean velocity profiles plotted in u/U M versus coordinates for U M = 12 m/s. 
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Figure 33. Boundary layer and internal layer thicknesses versus x for = 12 m/s. 
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Figure 36. Comparison of skin friction coefficient versus x for the rough-to-smooth, all- 
rough, and all-smooth surfaces for = 12 m/s. 




Figure 36. Mean velocity profiles for the rough-to-smooth surfaces compared with the 
smooth wall law of the wall for U M = 12 m/s. 
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Figure 37. Comparison of mean velocity profiles for the base-aligned and crest-aligned 
cases for U = 12 m/s in u/U M versus y/8 coordinates. 
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Figure 40. Comparison of skin friction coefficient versus x for the base-aligned and crest- 
aligned cases for = 58 m/s. 



CHAPTER V 

COMPARISON OF DATA AND PREDICTIONS 

The development of the turbulent boundary layer over a surface with a step 
change in surface roughness can be predicted with good accuracy using standard 
roughness and turbulence models. In this chapter, the discrete element method is 
described, and computations using this method are compared with selected data sets 
, from the preceding chapters. 

Discrete Element Method 

The two basic categories in which calculation efforts have fallen are (1) the 
equivalent sandgrain approach and (2) the discrete element approach. While both 
methods require experimental input, the equivalent sandgrain approach may require 
experimental data on the particular surface under consideration. On the other hand, 
the discrete element approach incorporates more basic physics of the process and uses 
a more generalized empirical input. It is therefore applicable to a broader spectrum 
of rough surfaces without requiring surface-specific experimental data. Since the 
discrete element approach is used for computations reported in this report, an 

overview of this method is presented next. 

The discrete element approach considers the mass, momentum and energy 
transport processes on the collection of individual roughness elements and the smooth 
surface between the elements. The basic idea is to formulate a system of partial 
differential equations that describes the mass, momentum and energy transport for 
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the flow over, around and between the roughness elements. In this method, the 
roughness effects are taken as an integral part of the flow problem and not (as with 
the equivalent sandgrain approach) as some ill-defined boundary condition. 

The discrete element method used in this work is formulated for roughness 
elements with three dimensional shapes (as opposed to transverse ribs) for which the 
element cross section can be approximated as circular at every height, y. Thus, the 
geometric description of the roughness element, d(y), is easily included in this predic- 
tion scheme. 

The steady (Reynolds-averaged), two-dimensional turbulent boundary layer 
equations presented here are for flow over a rough surface with roughness elements 
of uniform shape and spacing as derived by Taylor et al. (1984, 1985). The equations 


are: 

ipM) + ^-(p(V> = 0 

dx dy 


(8) 
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Empirical models for , -puV , the roughness element drag coefficient C D (y), and 

the roughness element Nusselt number Nu d (y) are necessary for closure. 

The blockage parameters (3 X and p y and the element shape descriptor d(y) 
require no empirical fluid mechanics input as they are determined solely from the 
geometry of the rough surface. Taylor et al. (1984) have shown for uniform three- 
dimensional roughness elements with circular cross-section that 


where for y > k, d(y) = 0 and both p x and p y become identically 1.0. 

The boundary conditions for the discrete element approach for rough wall flows 
are identical to those for smooth wall flows. The wall location (y = 0) is the smooth 
surface on which the roughness elements occur. At y = 0, u = v = 0 and H = H w . As 
y — » °°, u — » U„ and H -» H m . 

The numerical solution of the discrete element equations is obtained by finite 
difference solution of the transformed equations in the computational plane. The 
transformation, finite difference scheme, and program structure are described by 
Taylor et al. (1984). The steam wise derivative is approximated with a first-order 
backwards difference. The surface normal derivatives are replaced with second-order 
approximations which allow the spacing between grid points to vary with distance 
from the wall. This allows a concentration of nodes near the wall and below the crests 
of the elements. In this stretched grid the ratio of any two adjacent mesh lengths is 


a constant. 
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The solution is by an iterative marching, implicit method. The solution is 
known at station i and is sought at station i + 1. The implicit difference equations 
result in a tridiagonal coefficient matrix whose inverse is known and can be expressed 
algebraically (often referred to as the Thomas algorithm). Since the equations are 
non-linear, the system must be solved by iteration. A relaxation scheme is employed 

with a required residual < 0.01 percent. 

The solutions were obtained on finer and finer grids until differences were less 

than 1 percent in computed values of C f and St. In the transformed coordinates there 
were typically 120 grid points across the boundary layer with approximately 20-30 
grid points below the crests of the elements. The streamwise grid spacing was 

typically 1 cm. 

In addition, the codes were verified by comparisons with known solutions of 


smooth-wall laminar and turbulent flows. 

The "wall shear stress" is defined as the sum of the shear and drag forces on 
the wall in the mean flow direction divided by the plan area of the wall. The corre- 
sponding skin friction coefficient is then 



dy 


k 

i J_ f (p dC D u 2 )dy 
2 L 2 J 0 



( 12 ) 


and the Stanton Number is 
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St = 
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In order to solve Eqs. (8), (9) and (10), turbulence models for and^ 

and roughness models for C D and Nu d are required. Because of its proven predictive 
capability for attached boundary layer flows over smooth surfaces, the Prandtl mixing 
length model with van Driest damping and a constant turbulent Prandtl number is 
used for turbulence closure. Thus 


-pi/ V = pi 


3 u X 3 u 


3 u 


Oil I UUr j _ ^ 

1^-1 " 3 T 


3y ) 3y 




(14) 


where 

l m = 0.40y[l - exp(-y + /26)] ; l m < 0.095 
l m = 0.095; otherwise. 


(15) 


and 


-p iTA 7 = -1^ ^ 

Pr t 3y 


Pi 3 H 


(16) 


where Pr t = 0.9. Transition from laminar to turbulent flow is modeled using a 
streamwise intermittency factor 

Peff m ! _ e -3.61(x/x c -lf ( 17 ) 

Pt 


where x c is the specified point at which transition begins. 

After the surface changes abruptly from rough to smooth, the turbulence will 
not immediately take on the new equilibrium values for the smooth surface. The 
turbulent flow "remembers” the upstream history. The mixing length model presented 
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above does not account for this effect. Here, a fairly crude fix is adopted by using a 
lag model downstream of the step in roughness. This was done by taking a weighted 
average of the predicted local eddy viscosity, p*. with the eddy viscosity at some 


upstream location 


]i t = p i 0 (i-ri) + 


(18) 


where X = 55 0 with 5 0 being the boundary layer thickness at x 0 which was located 
0.02 m upstream of the interface between the surfaces. This idea was taken from a 
boundary layer shock wave interaction model of Shang, Hankey, and Law (1976). 

The roughness element C D and Nu d models are formulated as functions of the 
local element Reynolds number Re d = u(y)d(y)/v which includes roughness element 
size and shape information through d(y). As discussed in Taylor et al. (1984), the C D 
model which gave the best overall agreement with experimental data was 


logio(Cjr)) = - 0.125 logjo^e^) + 0.375 
This model has been tested for values of Re d up to 

25,000 (Taylor et al. (1984), Scaggs et al. (1988)) using many data sets. In particular, 
Scaggs et al. used eleven different rough surfaces, nine of which had uniform 
roughness elements and two of which were roughened nonuniformly. It was 
demonstrated that the roughness element drag coefficient model in the discrete 
element prediction approach gave excellent agreement with all of these data sets. 
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Consequently, this model was used unchanged for the predictions presented in this 
work. 

The roughness element energy transport model requires empirical input in the 
form of a Nusselt number, Nu d . Hosni, et al. (1989) developed the model 

Nu d = 1.7 Re°/ 9 Pr 0A (21) 

which is used in this work. This model has been tested up to Re d of about 2200. All 
of the experiments reported in this article have local element Reynolds numbers, Re d , 
which fall within the ranges of calibration listed above. 


Comparisons 

In the following comparisons all predictions and data correspond to the base 
aligned cases. 

Figures 41-46 show the comparisons between the data and the predictions for 
freestream velocities of 6, 12, 27, 43, 58, and 66 m/s. As stated before, the 6 m/s case 
in Figure 41 has a transitionally rough turbulent boundary layer over the rough 
portion of the surface. The other velocities have fully rough boundary layers over the 
rough portion. The discrete element method predictions are seen to be in good 
agreement for the 12 m/s and higher velocity cases (Figures 42-46). The rough portion 
Stanton number are somewhat underpredicted for the 6 m/s case, but the trends of 

this complicated case are reasonably predicted. 

Figure 47 shows the results of the velocity profile predictions for the 12 m/s 
case. These are the same conditions as the data given in Figure 24. Comparison of 
the two figures shows that the predictions agree well with the data. Figures 48 and 
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49 show the comparison for skin friction coefficients with freestream velocities of 12 
and 58 m/s respectively. The agreement is seen to be very good. 
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Figure 42. Comparison of base-aligned Stanton number predictions 
data for = 12 m/s. 
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Figure 48. Comparison of base-aligned skin friction coefficient predictions 
and data for = 12 m/s. 
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Figure 49. Comparison of base-aligned skin friction coefficient predictions 
and data for U = 58 m/s. 



CHAPTER VI 

SUMMARY AND CONCLUSIONS 


The purpose of this research effort was to investigate the effects of a step 
change in surface roughness from rough to smooth on heat transfer and skin friction 
in the turbulent boundary layer. Heat transfer and fluid mechanics data were 
collected for zero pressure gradient, incompressible, constant wall temperature air 
flows over the rough-to-smooth surface. These data were collected for both alignments 
with the smooth surface aligned with the base of the hemispherical roughness 
elements and the smooth surface aligned with the crests of the roughness elements. 
The rough-to-smooth data were compared to previously reported THTTF all-rough 
surface data (Hosni 1989) and all-smooth surface data (Coleman et al., 1988). 

For transitionally rough flow, the Stanton number after the rough-to-smooth 
interface quickly drops in a smooth, continuous fashion to the new smooth wall 
equilibrium value. For fully rough flows, on the other hand, the Stanton number 
undergoes an immediate drop to a value at or below the equivalent smooth-wall 
Stanton number at the same x-Reynolds number. The thermal boundary layer 
temperature profiles show the thermal boundary layer to slowly approach a smooth 
wall equilibrium profile. The alignment of the rough to smooth surfaces shows only 
a very small effect on the heat transfer between the base-aligned and crest aligned 
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cases. Any effect that is observed between the two alignment cases is limited to the 

region just downstream of the rough-to-smooth interface. 

The skin friction coefficient distributions behave very much like the Stanton 
number distributions. After the step change in surface roughness, the skin friction 
coefficient, C f , falls below the equivalent smooth-wall values before approaching the 
smooth-wall values. Mean velocity and turbulence intensity profiles show the flow to 
rapidly assume smooth-wall behavior in the near wall region, while requiring more 
distance to assume a complete smooth-wall behavior. Also seen from the fluid 
mechanics data is the emergence of two layers in the flow field downstream of the 
rough-to-smooth interface. In addition to the external boundary layer, an internal 
layer develops which separates the rough-wall flow characteristics from the rapidly 
developing smooth-wall flow. There are no discemable differences between the two 

alignment cases for the fluid mechanics data. 

Boundary layer computations using the standard discrete element method 

roughness model and mixing length turbulence model are shown to be in very good 
agreement with the data. 

These tests indicate that data obtained with smooth heat flux gauges on 
otherwise rough SSME fuel pump turbine blades could have large errors. We 
recommend that heat flux data taken in this fashion be used with caution. 
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APPENDIX A 


EXPERIMENTAL STANTON NUMBER DETERMINATION 
AND ITS UNCERTAINTY ANALYSIS 


Stanton Number Data Reduction Equation 

The Stanton number is the nondimensional convective heat transfer coefficient 

and may be defined as 


St = 


h 

P c p u„ 


( 22 ) 


where 

h is the convective heat transfer coefficient 
p is the density of the freestream air 
Cp is the specific heat of freestream air 
U m is the velocity of freestream air 

The convective heat transfer coefficient for heat transfer from a test plate to the air 
in the tunnel is defined as 


h = 


g 

A(T p - T a ) 


( 23 ) 


where 

q is the convective heat rate from the test plate 
A is the plate area 
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T p is the plate surface temperature 

T is the freestream air total temperature 

0 

Solving equation (9) for the convective heat transfer coefficient and substituting into 
equation (8) gives the Stanton number at each test plate as 


p C p U„A(T p - T a ) 

Since the convective heat transfer rate, q, will be determined, it is required 
that the corresponding radiation and conduction heat losses be known. Figure 41 
shows the modes of heat exchange from each plate. Application of an energy balance 

to a plate gives 


W = q + <7 C + Qr 


(25) 


where 

W is power supplied to the plate heater pad 

q r is the radiation heat loss rate 

n is the conduction heat loss rate 

Solving equation (11) for the convective heat transfer rate and substituting into 
equation (10) yields 


W - <? c - Qr 
p C p UJL(T p - TJ 

The radiation heat loss is modeled using 

q r = at A (T p 4 - T r 4 ) 

where 


(26) 


( 27 ) 


o 


is the Stefan-Boltzmann constant 


e 


is the test plate emissivity 

is the freestream recovery temperature 


93 


The conduction heat loss rate is modeled using 

q c , ( UA) pr (T p - T„a) * < U) PP (T P - T p .O <- (UA)p p (T p - T p .O < 28 ) 


where 


(UA) pr 


(UA) pp 


is the experimentally-determined effective conductance 
between a test plate and the side rails 
is the experimentally-determined effective conductance 
between two adjoining test plates 


T , is the side rail temperature at the axial location of the plate 

^rail 

rp ^ the temperature of the previous adjoining plate 

r p i the temperature of the next adjoining plate 

Substitution of equations (13) and (14) into equation (12) gives the final form 
of the data reduction equation for the experimentally determined Stanton number 


W - (UA) pr { Tp - Trail) - (UA)ppiiTp - Tp-t) + (^ P ~ r P+1 )) ~ q eA (T p 4 -7Y 4 ) 
St (29) 

Most of the variables involved in the experimentally determined Stanton 
number are shown explicitly in this expression. However, additional variable are used 
in the calculation of the static and total temperature of the freestream air and in the 
moist air property calculations for specific heat and density. The freestream air total 
and static temperatures are calculated using the measured recovery temperature and 
a recovery factor, R, (Eckert and Goldstein, 1976) for the probe: 
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T 0 = T r + (1-R) 


UJ 
2 c r 


( 30 ) 


T„~T r - (R) 


UJ 

2 C r 


( 31 ) 


The functional relationships for the moist air specific heat and density calculations are 


C p = CpiT^'T^^baro’Cpair’CpH^ 


( 32 ) 


p = p ( T „ , T w b , Pbaro ^ 


( 33 ) 


where 

T is the freestream air dry-bulb temperature 
T , is the freestream air wet-bulb temperature 
Pbaro is the barometric pressure 

q is the specific heat of dry air 
P air 

r is the specific heat of water vapor 
c ph 2 o 

In order to determine the Stanton number for each plate, the following sixteen 
variables must be measured or found from a reference source: 

• Plate heater power (W) 

• Recovery temperature (T r ) 

• Plate temperature (T p ) 

• Plate temperature (T^) 

• Plate temperature (T p+1 ) 

• Rail temperature (T rai] ) 
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• Wet-bulb temperature (T wb ) 

• Plate-to-rail effective conductance (UA) pr 

• piate-to-plate effective conductance (UA)pp 

• Plate area (A) 

• Barometric pressure (Pjjjro) 


• Specific heat of dry air C p air 


' ph 2 o 


• Specific heat of water vapor C p 

• Freestream air velocity (U M ) 

• Recovery factor (R) 

• Emissivity (e) 

Uncertainty Analysis Overview 

The detailed uncertainty analysis is consistent with the 1986 ANSI/ASME 
Standard on Measurement Uncertainty and follows the procedure of Coleman and 
Steele (1989). The uncertainty, U r , in an experimentally determined result, r, defines 
the interval r ± U r within which the true result lies with a 95 percent confidence. 
Here U r is the uncertainty in the result determined from a root-sum-square 
combination of the bias limit of the result, B r , and the precision limit of the result, P r 

U r = ( B r 2 + P r 2 ) 1/2 (34) 

For the Stanton number which is a function of sixteen variables, Xj 


r = r(Xi , Xo, ... ,X: ) 


(35) 


The propagation of precision limits, p Xj , of the measured variables and parameters 
into the result is given by 
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Pr- 


J 

E 


[±-Px, ] 


ax 


1/2 


(36) 


The propagation of bias limits, B Xj , of the measured variables and parameters into 

the result is given by 


11/2 


5r = 


J 

E 


(^B x ) 


ax,- 




(37) 


The primed terms on the right hand side of equation (23) represent the portions of 
bias limits of two variables which are perfectly correlated. For this case, the bulb 
thermistors used to measure various temperatures in the THTTF were all calibrated 
against the same standard. Therefore, a portion of the bias limit for each thermistor 
is perfectly correlated with a portion of the bias limits of the other thermistors. 

Previous work in the THTTF [Coleman et al. (1988)] has shown that most of 
the experiments in this facility are bias limit dominated. All of the P x , ’s are 

negligible when compared, to the bias limits. Therefore, 

P*-o. <38) 

An exception to this occurs at low freestream velocities (U M < 12 m/s), for which the 
heat transfer coefficients are low. At these conditions, the time constant of the 
THTTF is large enough that a tight steady state is hard to maintain. This is due 
mainly because of fluctuations in the line voltage to the test plate heater circuits and 
fluctuations in the temperature of the incoming cooling water. After observations of 
the St results for eight = 12 m/s replications with the smooth wall and three 
= 6 m/s replications with the rough wall, the 95% confidence estimate of P St = 3% 
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was determined. This is present because of system unsteadiness and not because of 
measurement uncertainty [Coleman et al. (1988)]. 

Applying equation (23) to equation (15) gives 

v . (J V * <#-> 2 *v , 2 * > 2 V 

Pi-1 


Pi.l 


, ( JiLfBr 2 . (_!®L) 2 B T * * (^) 2 Sw 2 * 

V dT~i Irail dT... I “ 6 dW dt/oo 


(39) 


dT^ 




, dSt ,2n 2 . c dSi In,..., 2 , <■ - ) 2 BrTTA't 2 ♦ (— ) 2 5c 2 

+ ( dC— ( dOMV } S(t/A) - W) pp (t7A) - * 


P Pi - 1 


P ^Pi+l 






* T r 


3S( w 3S( . n / Rt . / . 9( 3Si y 3S, _1B T 'B r ' 
^ "" "" 


Previous work in the THTTF has determined the following bias limits and 
nominal values for each of the sixteen variables involved in the calculation of the 
Stanton number. A detailed description of these estimates is given by Coleman et al. 
(1988) and Hosni et al. (1989). The values listed in parenthesis are the nominal 
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values for the 2.54 cm wide plates immediately downstream of the rough-to-smooth 
interface. 


Variable 

Bias Limit 

Nominal Values 

Plate Temperature 

0.14 C 

45 C 

Rail Temperature 

0.4 C 

45 C 

Recovery Temperature 

0.10 C 

30 C 

Wet-Bulb Temperature 

1.0 C 

27 C 

Barometric Pressure 

1.0 mm Hg 

760 mm Hg 

Recovery Factor 

0.09 

0.86 

Power 

0.9% 

20-150 W (3-18 W) 

Area 

0.03% 

464.5 cm 2 (116.1 cm 2 ) 

Freestream Velocity 

0.04% 

6-70 m/s 

C_ 

0.5% 

1.006 kJ/kg K 

Pair 

c_ 

0.5% 

1.86 kJ/kg K 

ph 2 o 

(UA) pr 

45% 

0.42 W/C (0.105 W/C) 

(UA) pp 

45% 

1.25 W/C 

Emissivity 

45% 

0.11 


Since the same calibration standard was used in the temperature calibration, 
the following bias limits are correlated 


= B t / = 0.04 C 

*■ rail 


A jitter program, as discussed by Coleman and Steele (1989), was used to 
determine the uncertainty estimate for each experimentally determined Stanton 
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number. The jitter program uses the Stanton number data reduction equation as a 
subroutine and approximates the partial derivatives in equation (25) using a finite 
difference scheme. Using this method, the Stanton number is first calculated using 
the original values of the sixteen variables for all twenty-four plates. Next, using a 
loop in the data reduction program, the partial derivative of Stanton number with 
respect to a particular variable is approximated by perturbing the variable a specified 
amount and recalculating the Stanton number for all of the plates. The partial 
derivative of Stanton with respect to the variable is the difference between the new 
Stanton number and the original Stanton number divided by the perturbation. The 
variables are then reset to their original values and the next variable is perturbed. 

After all of the partial derivatives have been calculated, the bias estimate for 
the uncertainty in the Stanton number is calculated directly from equation (25) using 
the partial derivatives and the previously stated bias limits. 

The overall uncertainty for the Stanton numbers reported in Chapter III 
ranged from 2 to 5 percent for the 10.2 cm wide plates, depending upon the flow 
conditions. For the 2.54 cm wide plate downstream of the rough-to-smooth interface, 
the uncertainty in the Stanton number was found to be from 3 to 11 percent, 
dependent upon the flow conditions. 
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Figure 50. Energy balance of a test plate. 



APPENDIX B 

TABULATION OF EXPERIMENTAL DATA 


Tabular listings of the experimental Stanton number and the skin friction 
coefficient data are presented in this Appendix. These data correspond to 
measurements of boundary layer flows over the rough-to-smooth test surface for both 
the base-aligned and crest-aligned cases. Stanton number data for both alignments 
are presented for freestream velocities of 6, 12, 27, 43, 58, and 66 m/s. Skin friction 
coefficient data for both alignments for freesteam velocities of 12 and 58 m/s are 

presented. 

As described in Appendix A, the Stanton number data were determined from 
an energy balance applied to each test plate. Plate 1 and plate 24 were considered as 
guard heaters and are not shown on the plots of heat transfer data. Their 
uncertainties are not reported in the listings. These data are marked with an in the 
following tabular listings. 

The first six digits of the "RUN" designation indicate the date in month, day, 
and year. The number after the dash is the run number for that particular day. 
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Stanton Number Data (Base-Aligned Case) 


Table 1. Stanton number data for the base-aligned case for LL = 6 m/s. 


RUN : 08179®- 1 


u - 

6.® 

m/s 

Pbar 

= 76.0 

Rho * 

1 . 17 

kg/n A 3 

Cp 

= 1.01 

Tr * 

Z7.G 

C 

T 

= 27. G 


PL 

x( n) 

Tu( C ) 

U(U> 

Qc( W ) 

1 

0.05 

44.8 

19. 1 

1 .8Z 

Z 

0. 15 

44.8 

12.1 

1.20 

3 

0.Z5 

44.7 

9.7 

0.54 

4 

0.36 

44.8 

1 1 .6 

0.44 

5 

0.4G 

44.8 

18.8 

0.40 

6 

0.5G 

44.7 

24.8 

-0. 1® 

7 

0.66 

44.6 

25.2 

-0. 1 1 

8 

0.76 

44.7 

23.7 

-0.01 

9 

0.86 

44.7 

22. 2 

-0.12 

1® 

0.93 

44.7 

4.5 

0.07 

1 1 

0.95 

44.7 

4.6 

-0.01 

1Z 

0.98 

44.7 

3.8 

0.01 

13 

1 .00 

44.7 

4. 1 

0. 10 

14 

1 .07 

44.6 

15.3 

0. 18 

15 

1 . 17 

44. 6 

15.5 

0.29 

16 

1 .27 

44.6 

14.4 

0.37 

17 

1.37 

44.8 

14.8 

0.61 

18 

1 .47 

44.7 

14.® 

0.35 

19 

1.57 

44.7 

14.0 

0.28 

Z0 

1 .68 

44.7 

13.7 

0.37 

Z1 

1.78 

44.8 

14.2 

0.84 

22 

1 .88 

44.8 

14.6 

0.84 

13 

1.98 

44.8 

13.7 

1 .04 

14 

2.08 

44.8 

25.5 

1 .74 


cn Hg 

Tub * 

18.9 C 

kJ/kg C 

Mu = 

1 . 86E-05 kg/m s 

C 

To =* 

27.6 C 


Qr(U) 

Rex 

St 

U <X> 

0.59 

1 . 90E+04 

0.00295 

0.0 * 

0.59 

5.71 E+04 

0.00182 

8.5 

0.59 

9.5ZE+04 

0.00153 

8.4 

0.59 

1 . 33E+05 

0.00188 

6.7 

0.59 

1 .71 E+05 

0.00315 

4.9 

0.59 

Z.09E+05 

0.00435 

4.3 

0.58 

Z.48E+05 

0.00441 

4.3 

0.59 

Z.86E+05 

0.00412 

4.4 

0.59 

3.24E+05 

0.00387 

4.5 

0.15 

3.48E+05 

0.00308 

10.7 

0. 15 

3.57E+05 

0.00319 

10.4 

0. 15 

3.G7E+05 

0.00256 

12.6 

0.15 

3.76E+05 

0.00275 

11.9 

0.58 

4.00E+05 

0.0026® 

5.4 

0.58 

4.38E+05 

0.00263 

5.6 

0.58 

4.76E+05 

0.00241 

5.8 

0.59 

5. 14E+05 

0.00242 

5.8 

0.59 

5.5ZE+05 

0.00233 

5.9 

0.59 

5.90E+05 

0.00234 

5.8 

0.59 

6.28E+05 

0.00227 

6.1 

0.59 

G.67E+05 

0.00226 

6.3 

0.59 

7.05E+05 

0.00233 

6.5 

0.59 

7.43E+05 

0.00213 

7.3 

0.59 

7.81E+05 

0.00410 

0.0 * 
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Table 2. Stanton number data for the base-aligned case for U„ - 12 m/s. 


RUN : 081590-1 

U = 1Z.0 n/s 

Rho = 1.17 kg/m A 3 

Tr - ZB. 5 C 


Pbar = 7G. 1 cm Hg 
Cp = 1.01 kJ/kg C 
T = ZB. 4 C 


Tub = 18.3 C 

Mu = 1.8BE-05 kg/m s 

To = ZB. 5 C 


PL 


1 

Z 

3 

4 

5 

6 

7 

8 
9 

10 

11 

1Z 

13 

14 

15 

16 

17 

18 
19 
Z0 
Z1 
ZZ 
Z3 
Z4 


x( n) 

Tu( C ) 

U<U) 

Qc(UI) 

Qr< U) 

Rex 

St 

U (X) 

0.05 

44. 1 

44.4 

0.00 

0.00 

0.00E+00 

0.00000 

0.0 * 

0. 15 

44.0 

7Z.9 

0.66 

0.60 

1.15E+05 

0.00614 

1.9 

0.25 

44. 1 

55. Z 

0.56 

0.60 

1.9ZE+05 

0.00461 

Z. 1 

0.36 

44. 1 

51 .4 

0.13 

0.60 

Z.69E+05 

0.00433 

Z. 1 

0.46 

44.0 

45.3 

-0.09 

0.60 

3.46E+05 

0.00383 

1.1 

0.56 

44.0 

4Z.B 

-0.Z5 

0.60 

4.Z3E+05 

0.00363 

2.3 

0.66 

44.0 

40. 1 

-0.Z0 

0.60 

5.00E+05 

0.00340 

Z.3 

0.76 

44.0 

38.8 

-0.48 

0.60 

5.77E+05 

0.00333 

Z . 4 

0.86 

44.0 

38. G 

-0. 17 

0.60 

6.54E+05 

0.003Z7 

2.4 

0.93 

44.0 

6.9 

-0.08 

0.15 

7.0ZE+05 

0.00Z34 

7.4 

0.95 

44.0 

6.5 

-0.13 

0. 15 

7.Z1E+05 

0.00ZZ4 

7.7 

0.98 

44. 1 

5.7 

-0.04 

0. 15 

7.40E+05 

0.00191 

8.7 

1 .00 

44. 1 

6. 1 

-0.05 

0. 15 

7.G0E+05 

0.00Z05 

8.2 

1 .07 

44. 1 

Z3.8 

0.07 

0.60 

8.08E+05 

0.00197 

4.4 

1 . 17 

44. 1 

Z4.5 

0.Z0 

0.60 

8.85E+05 

0.00Z03 

4.4 

1 .27 

44. 1 

Z3.Z 

0.Z8 

0.60 

9.61E+05 

0.00191 

4 . 5 

1.37 

44. 1 

Z3.4 

0.Z8 

0.60 

1 . 04E+06 

0.0019Z 

4.4 

1 .47 

44. 1 

11.1 

0. 1 1 

0.60 

1.1ZE+06 

0.00188 

4.5 

1 .57 

44. 1 

ZZ. 8 

0. 1 1 

0.60 

1 . 19E+0G 

0.00189 

4.4 

1 .68 

44. 1 

ZZ.6 

0.Z3 

0.60 

1 . Z7E+06 

0.00186 

4.5 

1.78 

44. 1 

11.1 

0.4Z 

0.60 

1 . 35E+06 

0.00185 

4. E 

1 .88 

44. 1 

Z3.6 

0.6Z 

0.60 

1.4ZE+06 

0.00191 

4 . 6 

1 .98 

44. 1 

11. S 

1 .04 

0.60 

1 . 50E+06 

0.00181 

4 . 9 

2.08 

44. 1 

36. 1 

0.00 

0.00 

0.00E+00 

0.00000 

0.0 * 
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Table 3. Stanton number data for the base-aligned case for U.. - 27 m/s. 


RUN : 

081G90-1 

U 

Z7.9 m/s 

Rho - 

1 . 17 kg/m A 3 

Tr * 

ZG.9 C 


Pbar * 7G. 1 cn Hg 
Cp = 1.01 kJ/kg C 
T = ZG.G C 


Tub - 18.3 C 

Mu = 1.8GE-05 kg/m s 

To - Z7.0 C 


PL 

x< n) 

Tu(C) 

U<U) 

Qc(W) 

Qr( U ) 

Rex 

St 

U (X) 

| 

0.05 

44.5 

1G6.G 

1.74 

0.60 

8.90E+04 

0.00610 

0.0 * 

z 

0. 15 

44.4 

1Z3.6 

1 .18 

0.60 

Z.67E+05 

0.00454 

1 .7 

3 

0.25 

44.4 

10Z.4 

0.60 

0.60 

4.45E+05 

0.00378 

1 .7 

4 

0.3G 

44.4 

99. Z 

0.4Z 

0.60 

6.Z3E+05 

0.00367 

1 .7 

5 

0.46 

44.4 

90. 1 

0.13 

0.60 

8.01E+05 

0.00334 

1 . 8 

6 

0.56 

44.5 

88. Z 

0.41 

0.60 

9.79E+05 

0.003Z4 

1 .8 

7 

0.66 

44.4 

83.8 

0.19 

0.60 

1 . 16E+06 

0.00309 

1.8 

8 

0.76 

44.4 

8Z. 1 

-0. 1 1 

0.60 

1 . 34E+06 

0.00305 

1.8 

g 

0.86 

44.5 

81 .9 

0.30 

0.60 

1 .51E+06 

0.00301 

1 .8 

10 

0.93 

44.4 

1 Z .7 

-0. 18 

0.15 

1.6ZE+06 

0.00190 

4. 1 

1 1 

0.95 

44.4 

1Z.0 

0.03 

0. 15 

1 . 67E+06 

0.00176 

4. 3 

1 2 

0.98 

44.5 

10.9 

0.07 

0.15 

1 .71 E+06 

0.00159 

4.7 

13 

1 .00 

44.5 

1 1.4 

0.01 

0. 15 

1 .76E+06 

0.00166 

4 . 5 

14 

1 .07 

44.5 

46. 1 

0.48 

0.60 

1 . 87E+06 

0.00167 

2.2 

15 

1 . 17 

44.5 

47. Z 

0.53 

0.60 

Z.05E+06 

0.00171 

Z . 2 

16 

1 .27 

44.5 

46. 1 

0.70 

0.60 

Z.Z3E+06 

0.00166 

1.1 

17 

1.37 

44.5 

45.9 

0.38 

0.60 

Z.40E+06 

0.00167 

2.2 

18 

1 .47 

44.5 

45.7 

0.5Z 

0.60 

Z.58E+06 

0.00165 

2.2 

19 

1 .57 

44.5 

45.7 

0.37 

0.60 

Z.76E+06 

0.00166 

1.1 

20 

1.68 

44.4 

45. 1 

0.35 

0.60 

Z.94E+06 

0.00164 

2.2 

2 1 

1 .78 

44.5 

46.0 

0.65 

0.60 

3.1ZE+06 

0.00166 

2.2 

22 

1 .88 

44. 5 

46.6 

0.83 

0.60 

3.Z9E+06 

0.00168 

2.3 

23 

1 .98 

44.5 

44.4 

1.03 

0.60 

3.47E+06 

0.00159 

2.4 

24 

2.08 

44.5 

6Z.3 

1.47 

0.60 

3.65E+06 

0.00ZZ4 

0.0 * 


105 


Table 4. Stanton number data for the base-aligned case for U M - 43 m/s. 


RUN : 081 690- Z 


u 

43.0 m/s 

Rho * 

1.16 kg/m*3 

Tr = 

Z8.9 C 


PL 

X< Fl) 

Tu(C) 

1 

0.05 

44.4 

Z 

0.15 

44. Z 

3 

0.Z5 

44.4 

4 

0.36 

44.4 

5 

0.46 

44.4 

6 

0.56 

44.3 

7 

0. GG 

44.4 

8 

0.76 

44.4 

9 

0.86 

44.4 

10 

0.93 

44.3 

1 1 

0.95 

44.3 

1Z 

0.98 

44.4 

13 

1.00 

44.3 

14 

1 .07 

44.3 

15 

1 . 17 

44.3 

16 

1 .11 

44.4 

17 

1.37 

44.4 

18 

1 .47 

44.4 

19 

1.57 

44.3 

Z0 

1 .68 

44.4 

Z1 

1.78 

44.4 

11 

1 .88 

44.4 

13 

1.98 

44. 3 

1 4 

Z .08 

44 . 4 


Pbar = 76.0 cm Hg 
Cp = 1.01 kJ/kg C 
T = Z8.1 C 


U(U) Qc( U ) Qr< U) 


Tub = 19.4 C 

Mu = 1.86E-05 kg/m s 

To = Z9.0 C 


Rex 


St 


U (X) 


Z16.7 

147.0 
1Z7.0 
1 ZB. 1 
1 14.6 
1 10.4 
107.3 

105.0 
105.0 

15.6 

14.4 

13.5 
13.8 

56.6 
58. 1 
57. Z 

58.1 

56.5 

56.7 

57. 1 

57.8 

57.6 

55.7 

11.1 


1.79 
0.66 
0.5Z 
0.Z7 
0.10 
-0.43 
- 0 . 1 1 
-0.34 
-0.ZZ 
- 0 . 1 1 
-0.18 
0. 10 
-0.Z6 
0. 16 
0.Z5 
0.36 
0.51 
0.18 
0.06 
0.49 
0.56 
0.68 
0.90 
1.63 


0.54 

0.53 

0.53 

0.54 

0.54 

0.53 

0.54 

0.54 

0.54 

0.13 

0.13 

0.13 

0.13 

0.53 

0.53 

0.53 

0.54 

0.53 

0.53 

0.54 

0.54 

0.53 

0.53 

0.53 


1 . 36E+05 

4.09E+05 

6.8ZE+05 

9.55E+05 

1 .Z3E+06 

1 . 50E+06 

1 . 77E+06 

Z.05E+06 

Z . 3ZE+06 

Z.49E+06 

Z . 56E+06 

Z.63E+06 

Z.69E+06 

Z.86E+06 

3. 14E+06 

3.41E+06 

3.68E+06 

3.96E+06 

4.Z3E+06 

4.50E+06 

4.77E+06 

5.05E+06 

5.3ZE+06 

5.59E+06 


0.00591 

0.00406 

0.00348 

0.00341 

0.00313 

0.00305 

0.00Z94 

0.00Z89 

0.00Z88 

0.0017Z 

0.00160 

0.00146 

0.00154 

0.00154 

0.00159 

0.00155 

0.00157 

0.00154 

0.00155 

0.00154 

0.00156 

0.00156 

0.00150 

0.00195 


0.0 • 
1.8 
1.8 
1.9 
1.9 
1.9 
1.9 
1 .9 

1.9 
3.7 

3.9 
4. 1 

4.0 
Z. 1 
Z. 1 
Z. 1 
Z. 1 

1.1 

1.1 

1 . 1 

Z. 1 

1.1 
1.1 
0.0 « 
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Table 5. Stanton number data for the base-aligned case for U„ - 58 m/s. 


RUN : 

081390-2 


U 

58.1 w/s 

Pbar =■ 76.Z 

Rho = 

1 .16 kg/n A 3 

Cp « 1.01 

Tr = 

31. Z C 

T = Z9.8 


PL x(n> Tu(C) U(U> Qc(U) 


1 0.05 44.3 ZZ4.5 1.81 

Z 0.15 44.3 158. Z 1.18 

3 0.Z5 44. Z 135.3 0.41 

4 0.3G 44. Z 133. Z 0.17 

5 0. 46 44. Z 1ZZ.3 -0.03 

G 0. 5G 44. 1 118.9 —0.37 

7 0. 66 44.1 115.1 -0.30 

8 0.76 44.1 113.0 -0.48 

9 0.86 44. Z 113.3 -0.37 

10 0.93 44. Z 16.7 -0.04 

11 0.95 44. Z 15.6 -0.12 

1 Z 0.98 44.2 14.3 -0.06 

13 1.00 44. Z 15.1 -0.08 

14 1.07 44. Z 60.9 -0.04 

15 1.17 44. Z 63.5 0.54 

1 6 1.Z7 44.1 61.1 0.2Z 

17 1.37 44.2 62. 4 0.5Z 

18 1.47 44.2 60.9 0.18 

19 1.57 44.3 61.8 0.43 

Z0 1.68 44. Z 60.8 0.34 

Z 1 1.78 44.1 61.1 0.34 

Z2 1.88 44.2 61.9 0.65 

23 1.98 44. Z 60.2 0.97 

24 Z . 08 44. Z 74.7 1.48 


cm Hg 

Tub 

■> Z0.0 C 

kJ/kg C 

Mu 

= 1.86E-05 kg/n s 

C 

To 

= 31.5 C 


Qr< U) 

Rex 

St 

U <X) 

0.46 

1 . 84E+05 

0.00546 

0.0 * 

0.45 

5.51E+05 

0.00386 

Z . 4 

0.45 

9. 18E+05 

0.00333 

Z . 4 

0.45 

1.Z8E+06 

0.003Z9 

Z . 5 

0.45 

1.65E+06 

0. 00303 

Z . 5 

0.45 

Z.0ZE+06 

0.00Z96 

Z.6 

0.45 

Z . 39E+06 

0.00Z86 

Z.5 

0.45 

Z.75E+06 

0.00Z8Z 

Z.6 

0.45 

3.1ZE+06 

0.00Z81 

Z.6 

0.11 

3.35E+06 

0.00165 

3.9 

0.11 

3.44E+06 

0.00155 

4. 1 

0.11 

3.53E+0G 

0.0014Z 

4.3 

0.11 

3 . 63E+06 

0.00150 

4. 1 

0.45 

3.85E+06 

0.001 50 

1.1 

0.45 

4.ZZE+06 

0.00155 

Z.6 

0.45 

4.59E+06 

0.00151 

1.1 

0.45 

4.96E+06 

0.0015Z 

Z.6 

0.45 

5.3ZE+06 

0.00150 

Z.6 

0.45 

5.69E+06 

0.00150 

Z.6 

0.45 

6.06E+06 

0.00149 

2.6 

0.45 

6.4ZE+06 

0.00150 

Z.6 

0.45 

6.79E+06 

0.00151 

Z.6 

0.45 

7.16E+06 

0.00146 

Z.6 

0.45 

7.53E+06 

0.00181 

0.0 * 
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Table 6. Stanton number data for the base-aligned case for = 66 m/s. 


RUN : 
U - 
Rho ■ 
Tr * 

081490-1 
66 .3 n/a 
1.15 kg/n'' 3 
33.5 C 

Pbar 

Cp 

T 

= 76. 1 
= 1.01 
= 31.6 

cn Hg 
kJ/kg C 
C 

Tub 

Mu 

To 

PL 

x< n ) Tu< C ) 

U<U> 

Qc( U) 

Qr< U) 

Rex 


1 

0.05 

44.3 

205.7 

2 

0.15 

44.2 

145.5 

3 

0.25 

44.2 

125. 1 

4 

0.36 

44.2 

123. 1 

5 

0.46 

44.2 

113.7 

6 

0.56 

44.2 

1 10.5 

7 

0.66 

44.2 

106.6 

8 

0.76 

44.2 

104.5 

9 

0.86 

44.3 

105.3 

10 

0.93 

44.2 

15.2 

11 

0.95 

44.2 

14.3 

12 

0.98 

44.2 

13.0 

13 

1 .00 

44.2 

14.0 

14 

1 .07 

44.2 

56.6 

15 

1 . 17 

44.3 

58.6 

16 

1.27 

44.3 

57.5 

17 

1.37 

44.2 

57.0 

18 

1 .47 

44.3 

56.8 

19 

1 .57 

44.3 

56.4 

20 

1.68 

44.2 

55.6 

21 

1 .78 

44.2 

56.6 

22 

1.88 

44.2 

56.7 

23 

1.98 

44.2 

55.0 

24 

2.08 

44.2 

66.8 


1.98 

0.38 

Z.08E+05 

1.32 

0.38 

6.Z4E+05 

0.58 

0.38 

1 .04E+06 

0.22 

0.38 

1 .4GE+06 

0.12 

0.38 

1 . 87E+06 

-0.30 

0.38 

Z.Z9E+06 

-0.22 

0.38 

Z.71E+06 

-0.41 

0.38 

3.1ZE+06 

-0.18 

0.38 

3.54E+0G 

-0. 1 1 

0.09 

3.80E+06 

-0.07 

0.09 

3.90E+0G 

-0. 14 

0.09 

4.01E+06 

-0. 10 

0.09 

4. 1 1E+06 

0.16 

0.38 

4.37E+06 

0.54 

0.38 

4.79E+06 

0.75 

0.38 

5.20E+06 

0.36 

0.38 

5.GZE+06 

0.57 

0.38 

6.03E+06 

0.34 

0.38 

6.45E+06 

0.26 

0.38 

6.87E+06 

0.53 

0.38 

7.Z8E+06 

0.50 

0.38 

7.70E+06 

0.90 

0.38 

8.12E+06 

1.38 

0.38 

8.53E+06 


21.1 C 

1.8GE-0S kg/n s 
33.8 C 


St U (X) 


0.00538 

0.0 * 

0.00382 

3.3 

0.00330 

3.4 

0.00326 

3.4 

0. 00300 

3.4 

0.00294 

3.5 

0.00283 

3.5 

0.00Z78 

3.5 

0.00278 

3.5 

0.00162 

4.8 

0.00152 

4.9 

0.00139 

5.2 

0.00149 

5.0 

0.00149 

3.5 

0.00153 

3.4 

0.00149 

3.4 

0.00150 

3.5 

0.00147 

3.4 

0.00147 

3.5 

0.00146 

3.5 

0.00148 

3.4 

0.00149 

3.5 

0.00143 

3.4 

0.00174 

0.0 * 
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Stanton Number Data (Crest-Aligned Case) 


Table 7. Stanton number data for the crest-aligned case for = 6 m/s. 


RUN : 102390-1 


u 

6.0 

m/s 

Pbar 

= 75.9 

cm Hg 

Tub * 

16.7 C 

Rho * 

1 .18 

kg/m* 3 

Cp 

= 1.01 

kJ/kg C 

Mu 

1.86E-05 kg/m s 

Tr = 

25. Z 

C 

T 

= 25. Z 

C 

To 3 

25.2 C 


PL 

x< n) 

Tu( C ) 

U<U) 

Qc( W) 

Qr(U) 

Rex 

St 

U <X> 

1 

0.05 

44. 1 

19.4 

1.31 

0.64 

1 .9ZE+04 

0.00279 

0.0 * 

2 

0. 15 

44. 1 

12.3 

0.73 

0.64 

5.75E+04 

0.00175 

7.2 

3 

0.25 

44. 1 

10.2 

0.29 

0.64 

9.58E+04 

0.00149 

7.4 

4 

0.36 

44.2 

1Z.6 

0. 19 

0.64 

1 . 34E+05 

0.00187 

6.2 

5 

0.46 

44. 1 

20.4 

-0.27 

0.64 

1 . 7ZE+05 

0.00321 

4.6 

6 

0.56 

44. Z 

28. 1 

0.02 

0.65 

Z. 1 1E+05 

0.00436 

4. 1 

7 

0.66 

44.2 

27.3 

-0.22 

0.64 

2 . 49E+05 

0.00429 

4.2 

8 

0.76 

44. 1 

Z5.3 

-0.38 

0.64 

Z . 87E+05 

0.00400 

4.3 

9 

0.86 

44.2 

23.7 

-0.15 

0.64 

3.Z6E+05 

0.00371 

4.4 

10 

0.93 

44. 1 

5. 1 

-0.27 

0.16 

3.50E+05 

0.00333 

9.4 

1 1 

0.95 

44.2 

4.9 

0.21 

0.16 

3.59E+05 

0.00292 

10.5 

12 

0.98 

44. 1 

4. 1 

-0.28 

0.16 

3.69E+05 

0.00271 

1 1.2 

13 

1 .00 

44.2 

4.4 

0.06 

0. 16 

3.78E+05 

0.00268 

11.1 

14 

1 .07 

44. 1 

16.7 

-0.09 

0.64 

4.0ZE+05 

0.00258 

5. 1 

15 

1 . 17 

44.2 

16.8 

0.37 

0.64 

4.41E+05 

0.00253 

5.2 

16 

1.27 

44. 1 

15.7 

0. 16 

0.64 

4.79E+05 

0.00238 

5.4 

17 

1.37 

44. 1 

15.6 

0.04 

0.64 

5.17E+05 

0.00239 

5.4 

18 

1 .47 

44. 1 

15.3 

0 . 11 

0.64 

5.5SE+05 

0.00231 

5.5 

19 

1 .57 

44.0 

15.0 

-0. 10 

0.64 

5.94E+05 

0.00232 

5.4 

20 

1 .68 

44.0 

14.7 

0.06 

0.64 

6.3ZE+05 

0.00225 

5.5 

21 

1 .78 

44. 1 

15.0 

0.36 

0.64 

6.70E+05 

0.00225 

5.6 

22 

1.88 

44. 1 

15.6 

0.45 

0.64 

7.09E+05 

0.00233 

5.6 

23 

1 .98 

44. 1 

15.0 

0.80 

0.64 

7.47E+05 

0.00217 

6.2 

24 

2.08 

44. 1 

21.8 

1 .18 

0.64 

7.85E+05 

0.00321 

0.0 * 
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Table 8. Stanton number data for the crest-aligned case for - 12 m/s. 


RUN : 

101Z90-1 



U 

1Z.Z 

m/s 

Pbar 

= 75.7 

Rho = 

1 . 17 

kg/m A 3 

Cp 

= 1.01 

Tr = 

Z5.7 

C 

T 

= 25. G 


PL 

x< m> 

Tu(C> 

U(U> 

Qc( U ) 

1 

0.05 

44. 1 

48.5 

0.00 

Z 

0. 15 

43.9 

75.2 

0.79 

3 

0.Z5 

43.9 

56.2 

0.37 

4 

0.3G 

43.9 

52.9 

0.42 

5 

0.4G 

43.9 

47.3 

0.12 

6 

0.5G 

43.9 

44.3 

-0.18 

7 

0.G6 

44.0 

42.1 

0. 1 1 

8 

0.76 

43.9 

40.6 

-0.08 

9 

0.86 

43.9 

39.3 

-0.26 

10 

0.93 

43.9 

7.5 

-0.02 

1 1 

0.95 

43.9 

6.8 

0. 10 

12 

0.98 

43.9 

6.0 

-0.20 

13 

1 .00 

43.9 

6.3 

0.09 

14 

1 .07 

43.9 

24.9 

0.06 

15 

1 . 17 

43.9 

25.5 

0.21 

1 G 

1 .27 

44.0 

24.4 

0.44 

17 

1.37 

44.0 

24.3 

0. 10 

18 

1 .47 

44.0 

23.9 

0.21 

19 

1.57 

44.0 

24.2 

-0.01 

20 

1 .68 

44.0 

23.6 

0. 12 

Z1 

1 .78 

44.0 

24. 1 

0.36 

22 

1.88 

44.0 

24.7 

0.48 

13 

1.98 

43.9 

23.3 

0.49 

24 

2.08 

43.9 

31 .0 

0.00 


cm Hg 

Tub = 

16.1 C 

kJ/kg C 

Mu = 

1.86E-05 kg/ pi s 

C 

To = 

25.7 C 


Qr( W) 

Rex 

St 

U (X) 

0.00 

0. 00E+00 

0.00000 

0.0 * 

0.62 

1.17E+05 

0.00602 

1.9 

0.6Z 

1 . 95E+05 

0.00451 

2.0 

0.62 

Z.73E+05 

0.00422 

Z. 1 

0.62 

3.51E+05 

0.00379 

2.1 

0.62 

4.Z9E+05 

0.00358 

1.1 

0.62 

5.07E+05 

0.00336 

2.3 

0.62 

5.85E+05 

0.00327 

2.3 

0.62 

6.63E+05 

0.00318 

2.3 

0. 15 

7.1ZE+05 

0.00240 

7.0 

0. 16 

7.3ZE+05 

0.00212 

7.7 

0.15 

7.51 E+05 

0.00197 

8.2 

0.15 

7.71E+05 

0.00198 

8. 1 

0.62 

8.19E+05 

0.00197 

4.3 

0.6Z 

8.97E+05 

0.00201 

4.3 

0.62 

9.75E+05 

0.00189 

4.4 

0.62 

1 . 05E+06 

0.00192 

4.3 

0.62 

1 . 1 3E+06 

0.00187 

4.4 

0.62 

1 .21E+06 

0.00191 

4.3 

0.62 

1 .Z9E+06 

0.00185 

4.4 

0.62 

1 . 37E+06 

0.00188 

4.4 

0.62 

1 . 44E+06 

0.00192 

4.4 

0.62 

1 . 5ZE+06 

0.00182 

4.6 

0.00 

0.00E+00 

0.00000 

0.0 * 
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Table 9. Stanton number data for the crest-aligned case for U.. - 27 m/s. 


RUN : 

101790-1 

U “ 

Z7.8 

m/s 

Rho = 

1 . 17 

kg/m* 3 

Tr ■ 

ZG.4 

C 


Pbar = 75.9 cm Hg 
Cp * 1.01 kJ/kg C 
T = ZG.1 C 


Tub - 18.9 C 

Mu = 1.8GE-05 kg/m s 

To * ZG.5 C 


PL 

x( m) 

Tu< C ) 

U(U) 

Qc( U) 

1 

0.05 

44. 1 

164. 1 

1.58 

Z 

0. 15 

44.0 

1ZZ.5 

0.86 

3 

0.Z5 

44.0 

101.8 

0.4Z 

4 

0. 36 

44.0 

96. 3 

0.Z5 

5 

0.4G 

44.0 

90.3 

0.15 

G 

0.56 

44.0 

8G.4 

-0.16 

7 

0.GG 

44.0 

8Z.8 

0.0Z 

8 

0.7G 

44.0 

80. 1 

-0.38 

9 

0.8G 

44.0 

79.9 

-0.04 

10 

0.93 

44.0 

13.3 

-0. 14 

11 

0.95 

44. 1 

1 1.8 

0.08 

1Z 

0.98 

44. 1 

10.9 

-0.07 

13 

1 .00 

44. 1 

1 l.Z 

0. 1Z 

14 

1.07 

43.9 

44.9 

-0.ZZ 

15 

1.17 

44.0 

46.4 

0.39 

1 G 

1 .11 

44.0 

45.0 

0.33 

17 

1.37 

43.9 

44.9 

-0.08 

18 

1 .47 

44.0 

44.5 

0. 16 

19 

1.57 

43.9 

44.8 

-0.06 

Z0 

1 .68 

43.9 

44.3 

0.09 

Z1 

1.78 

44.0 

44.7 

0.Z4 

ZZ 

1 .88 

43.9 

45.3 

0.43 

13 

1.98 

43.9 

43.9 

0.56 

Z4 

Z.08 

44.0 

53.0 

1.39 


Qr(W) 

Rex 

St 

U (X) 

0.60 

8.86E+04 

0.00600 

0.0 * 

0.60 

Z.66E+05 

0.00451 

1 .7 

0.60 

4.43E+05 

0.00376 

1.7 

0.60 

6.Z0E+05 

0.00355 

1 .7 

0.60 

7.98E+05 

0.00333 

1.8 

0.60 

9.75E+05 

0.003Z0 

1.8 

0.60 

1 . 15E+06 

0.00306 

1.8 

0.60 

1 . 33E+06 

0.00Z99 

1.8 

0.60 

1 .51E+0G 

0.00Z95 

1.8 

0.15 

1.6ZE+06 

0.00198 

4.0 

0.15 

1 . 66E+06 

0.0017Z 

4.4 

0. 15 

1.71E+06 

0.00161 

4.6 

0. 15 

1.75E+06 

0.00163 

4.6 

0.59 

1.86E+06 

0.00167 

1.1 

0.60 

Z.04E+06 

0.00169 

1.1 

0.60 

Z.ZZE+06 

0.00164 

1.1 

0.59 

Z.39E+0G 

0.00166 

1.1 

0.60 

Z.57E+06 

0.00163 

1.1 

0.60 

Z.75E+0G 

0.00165 

1.1 

0.60 

Z.9ZE+06 

0.00163 

1.1 

0.60 

3. 10E+06 

0.00164 

1.1 

0.60 

3.Z8E+06 

0.00166 

1.1 

0.60 

3.46E+06 

0.00160 

1.3 

0.60 

3.63E+06 

0.00190 

0.0 * 
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Table 10. St an ton number data for the crest-aligned case for U.. - 43 m/s. 


RUN : 101 89®- 1 

U 3 42.9 n/a 

Rho = 1 . 17 kg/n A 3 

Tr = 28. B C 


PL 

x< n) 

Tu( C ) 

1 

0.05 

44.2 

2 

0. 15 

44.4 

3 

0.25 

44.4 

4 

0.36 

44.3 

5 

0.46 

44.3 

6 

0. 5G 

44.4 

7 

0.66 

44.4 

8 

0.76 

44.4 

9 

0.86 

44.4 

1® 

0.93 

44.3 

1 1 

0.95 

44.4 

12 

0.98 

44.4 

13 

1 .00 

44.4 

14 

1 .07 

44.3 

15 

1.17 

44.4 

16 

1 .27 

44.4 

17 

1.37 

44 . 3 

18 

1 .47 

44.4 

19 

1.57 

44.4 

20 

1 .68 

44.3 

21 

1.78 

44.4 

22 

1.88 

44.4 

23 

1.98 

44.4 

24 

2.08 

44.4 


Pbar = 76. 3 
Cp - 1.01 
T =* 27.8 


U<U) 

Qc(UI) 

219.3 

1.30 

151.6 

1 .28 

129.0 

0.62 

120.8 

0. 14 

1 15.8 

0.03 

113.3 

0.15 

108.4 

-0.14 

106.4 

-0.03 

104.5 

-0. 13 

16.9 

-0.13 

14.9 

0. 1® 

13.8 

-0.07 

14.3 

0.13 

57.8 

-0.18 

60.6 

0.68 

58.2 

0.33 

58.2 

-0.05 

58. 1 

0.34 

58.3 

0.12 

56.9 

-0.05 

58.5 

0.47 

58.7 

0.36 

57.9 

0.93 

65.6 

1 . 17 


cn Hg Tub * 

kJ/kg C Mu = 

C To = 


Qr( U) Rex 


0.54 1.37E+05 
0.54 4.10E+05 
0.54 G.84E+05 
0.54 9.58E+05 
0.54 1.23E+0G 
0.54 1.50E+0G 
0.54 1.78E+06 
0.54 2.05E+0G 
0.54 2.33E+0G 
0. 14 2.50E+0G 
0.14 2.57E+0G 
0.14 Z.63E+06 
0.14 2.70E+0G 
0.54 2.87E+0G 
0.55 3.15E+0G 
0.54 3.42E+0G 
0.54 3.G9E+0G 
0.54 3.97E+0G 
0.54 4.24E+0G 
0.54 4.51E+06 
0.54 4.79E+06 
0.54 5.0GE+0G 
0.54 5.34E+0G 
0.54 5.G1E+0G 


1 G.7 C 

1.8GE-05 kg/n s 
28.7 C 


St 

U (X) 

0.00596 

0.0 * 

0.00407 

1.8 

0.00347 

1.8 

0.00328 

1.8 

0.00314 

1.8 

0.00305 

1.8 

0.00293 

1.9 

0.00287 

1.9 

0.00283 

1.9 

0.00183 

3.5 

0.00159 

3.8 

0.00149 

3.9 

0.00152 

3.9 

0.00157 

2.1 

0.0016® 

2.1 

0.00155 

Z. 1 

0.00157 

2.1 

0.00155 

2.1 

0.00156 

2.1 

0.00153 

2.1 

0.00156 

2.1 

0.00157 

2.1 

0.00153 

2.1 

0.00174 

0.0 * 
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Table 11. Stanton number data for the crest-aligned case for U.. = 58 m/s. 


RUN : 101 69©- 1 


u « 

58.0 

n/s 

Pbar 

“ 76.0 

cm Hg 

Tub * 

18.9 C 

Rho = 

1 .16 

kg/n A 3 

Cp 

- 1.01 

kJ/kg C 

Mu » 

1.86E-05 kg/n s 

Tr * 

31 . 1 

C 

T 

= Z9.7 

C 

To 33 

31.3 C 


PL 

x ( m) 

Tu(C) 

U(U) 

Qc< U) 

Qr(W) 

Rex 

St 

U (X) 

1 

0.05 

44.4 

Z30.0 

Z.Z6 

0.47 

1 . 83E+05 

0.00548 

0.0 * 

Z 

0. 15 

44.5 

16Z.9 

1.76 

0.47 

5.50E+05 

0.00386 

Z.3 

3 

0.Z5 

44.4 

139.Z 

1.01 

0.47 

9. 16E+05 

0.0033Z 

Z . 3 

4 

0.36 

44.5 

134.0 

0.93 

0.47 

1 . Z8E+06 

0.00318 

Z.3 

5 

0.46 

44.4 

1Z6.Z 

0.48 

0.47 

1 . 65E+06 

0.0030Z 

Z.3 

6 

0.56 

44.4 

1Z3.Z 

0.31 

0.46 

Z.01E+06 

0.00Z96 

Z . 4 

7 

0.66 

44.4 

118.8 

0.Z4 

0.46 

Z . 38E+06 

0.00Z86 

Z.4 

8 

0.76 

44.4 

1 16.0 

0. 14 

0.46 

Z.75E+06 

0.00Z79 

Z . 4 

9 

0.86 

44.4 

1 15. 1 

0.40 

0.47 

3.11E+06 

0.00Z75 

Z.3 

10 

0.93 

44.3 

18.0 

-0. 18 

0.1Z 

3. 34E+06 

0.00175 

3.8 

1 1 

0.95 

44.4 

16.0 

0.19 

0. 1Z 

3.43E+06 

0.00151 

3.9 

1Z 

0.98 

44.4 

14.9 

-0.07 

0.1Z 

3.53E+06 

0.00144 

4. 1 

13 

1.00 

44.4 

15.5 

0.Z3 

0. 1Z 

3.6ZE+06 

0.00147 

3.9 

14 

1 .07 

44.3 

63.0 

0.15 

0.46 

3.85E+06 

0.0015Z 

Z . 6 

15 

1 . 17 

44.4 

65.5 

0.69 

0.46 

4.Z1E+06 

0.00156 

Z.5 

16 

1 .Z7 

44.4 

63.7 

0.63 

0.46 

4.58E+06 

0.00151 

Z.5 

17 

1.37 

44.4 

64. Z 

0.61 

0.46 

4.95E+06 

0.0015Z 

Z.S 

18 

1.47 

44.4 

6Z.3 

0.38 

0.46 

5.31E+06 

0.00149 

Z.5 

19 

1.57 

44.4 

63.6 

0.40 

0.46 

5.68E+06 

0.00151 

Z.5 

Z0 

1 .68 

44.4 

6Z.4 

0.30 

0.46 

6.04E+06 

0.00149 

Z.5 

Z1 

1 .78 

44.4 

63. 1 

0.54 

0.46 

6.41E+06 

0.00150 

Z.5 

ZZ 

1.88 

44.4 

64. Z 

0.78 

0.46 

6.78E+06 

0.0015Z 

Z.5 

Z3 

1.98 

44.3 

61 .7 

0.7Z 

0.46 

7. 14E+06 

0.00147 

Z.5 

Z4 

Z . 08 

44.4 

70.6 

1 .6Z 

0.47 

7.51E+06 

0.00165 

0.0 * 
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Table 12. Stanton number data for the crest-aligned case for = 66 m/s. 


RUN : 10ZZ90-1 

U * 66. Z m/s 

Rho - 1.15 kg/ m A 3 

Tr = 3Z.7 C 


PL 

x( m) 

Tu(C) 

1 

0.05 

44. 1 

Z 

0. 15 

44. 1 

3 

0.Z5 

44. Z 

4 

0.36 

44. Z 

5 

0.46 

44. 1 

6 

0.56 

44. 1 

7 

0.66 

44. 1 

8 

0.76 

44. 1 

9 

0.86 

44. Z 

10 

0.93 

44. 1 

11 

0.95 

44. Z 

1Z 

0.98 

44. 1 

13 

1 .00 

44. Z 

14 

1 .07 

44. 1 

15 

1 . 17 

44. Z 

16 

1 . Z7 

44. Z 

17 

1.37 

44. 1 

18 

1 .47 

44. Z 

19 

1.57 

44. 1 

Z0 

1.68 

44. 1 

Z1 

1 .78 

44. 1 

ZZ 

1.88 

44. 1 

Z3 

1.98 

44. 1 

Z4 

Z.08 

44. Z 


Pbar = 75.8 
Cp =1.01 
T = 30.8 


U(U) Qc( W ) 


Z 14. 9 1.65 

154.4 0.93 

133.6 0.65 

1Z5.5 0.31 

1Z0.Z -0.07 

117.5 -0.06 

113.1 -0.11 

111.1 - 0.11 

109.8 0.08 

17.3 -0.18 

15. Z 0.1Z 

14. Z -0.13 

14.9 0.19 

60.7 0.09 

6Z.8 0.68 

61. Z 0.58 

60.8 0.Z4 

59.9 0.40 

60.1 0.13 

59. Z 0.18 

59.8 0.3Z 

61.0 0.41 

59. Z 0.58 

66.0 1.10 


cm Hg Tub = 
kJ/kg C Mu = 
C To = 


Qr< U) Rex 


0.40 Z.08E+05 
0.40 6.Z3E+05 
0.41 1 . 04E+06 
0.41 1.45E+06 
0.40 1.87E+06 
0.40 Z . Z8E+06 
0.40 Z.70E+06 
0.40 3. 1 1E-l-06 
0.40 3. 53E+06 
0.10 3.79E+06 
0.10 3. 89E+06 
0.10 4.00E+06 
0.10 4.10E+06 
0.40 4. 36E+06 
0.41 4.78E+06 
0.41 5.19E+06 
0.40 5.61E+06 
0.40 6.0ZE+06 
0.40 6.44E+06 
0.40 6. 85E+06 
0.40 7.Z7E+06 
0.40 7.68E+06 
0.40 8. 10E+06 
0.40 8.51E+06 


Z0.0 C 

1.86E-05 kg/m s 
33.0 C 


St U (X) 


0.0053Z 0.0 * 

0.00383 3.0 

0.003Z8 3.0 

0.00309 3.1 

0.00Z99 3.Z 

0.00Z9Z 3.1 

0.00Z8Z 3. Z 

0.00Z76 3.Z 

0.00Z7Z 3.1 

0.00174 4.4 

0.00149 4.5 

0.0014Z 4.8 

0.00145 4.5 

0.00150 3.3 

0.00153 3.1 

0.00149 3.1 

0.00150 3. Z 

0.00147 3. Z 

0.00149 3. Z 

0.00147 3.Z 

0.00147 3. Z 

0.00150 3. Z 

0.00145 3. Z 

0.00160 0.0 * 
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Skin Friction Coefficient Data (Base-Aligned Case) 

Table 13. Skin friction coefficient data for the base-aligned case for U,. = 12 m/s. 

The Cf value at x = 0.86 m was obtained using the hot-wire technique 
(Hosni, 1989), all other values were measured using the Preston tube 
technique. 


x (m) Cf 


0.86 

0.00616 

0.97 

0.00257 

1.07 

0.00292 

1.17 

0.00290 

1.27 

0.00293 

1.37 

0.00300 

1.58 

0.00302 

1.78 

0.00303 

1.98 

0.00301 
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Table 14. 


Skin friction coefficient data for the base-aligned case for U.. = 58 m/s. 
The Cf value at x = 0.86 m was obtained using the hot-wire technique 
(Hosni, 1989), all other values were measured vising the Preston tube 
technique. 


x (m) 


C f 


0.86 

0.00577 

0.97 

0.00188 

1.07 

0.00210 

1.17 

0.00222 

1.27 

0.00230 

1.37 

0.00234 

1.58 

0.00237 

1.78 

0.00239 

1.98 

0.00237 
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Skin Friction Coefficient Data (Crest- Aliened Case) 

Table 15. Skin friction coefficient data for the crest-aligned case for U„ = 12 m/s. 

The C f value at x = 0.86 m was obtained using the hot-wire technique 
(Hosni, 1989), all other values were measured using the Preston tube 
technique. 


x (m) C f 


0.86 

0.00616 

0.97 

0.00250 

1.07 

0.00279 

1.17 

0.00289 

1.27 

0.00295 

1.37 

0.00296 

1.58 

0.00299 

1.78 

0.00303 

1.98 

0.00303 
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Table 16. 


S irin fric tio n coefficient data for the crest-aligned case for U„ — 58 m/s. 
The C f value at x = 0.86 m was obtained using the hot-wire technique 
(Hosni, 1989), all other values were measured using the Preston tube 
technique. 


X (m) 

c f 

0.86 

0.00577 

0.97 

0.00191 

1.07 

0.00216 

1.17 

0.00224 

1.27 

0.00231 

1.37 

0.00234 

1.58 

0.00237 

1.78 

0.00239 

1.98 

0.00240 


